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We start by reviewing the existing literature on the creation of black holes during high- 
energy particle collisions, both in the absence and in the presence of extra, compact, 
spacelike dimensions. Then, we discuss in detail the properties of the produced higher- 
dimensional black holes, namely the horizon radius, temperature and life-time, as well as 
the physics that governs the evaporation of these objects, through the emission of Hawking 
radiation. We first study the emission of visible Hawking radiation on the brane: we derive 
a master equation for the propagation of fields with arbitrary spin in the induced-on-the- 
brane black hole background, and we review all existing results in the literature for the 
emission of scalars, fermions and gauge bosons during the spin-down and Schwarzschild 
phases of the life of the black hole. Both analytical and numerical results for the greybody 
factors and radiation spectra are reviewed and exact results for the number and type of 
fields emitted on the brane as a function of the dimensionality of spacetime are discussed. 
We finally study the emission of Hawking radiation in the bulk: greybody factors and 
radiation spectra are presented for the emission of scalar modes, and the ratio of the 
missing energy over the visible one is calculated for different values of the number of extra 
dimensions. 
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1. Introduction 

A few years ago, in 1998, a novel idea was proposed according to which the so- 
called hierarchy problem - in other words, our difficulty in answering the question 
of why the characteristic scale of gravity, Mp ~ 10 19 GeV, is 16 orders of magnitude 
larger than the Electro- Weak scale, Mew ~ 1 TeV - could be solved by assuming 
the existence of extra dimensions in the Universe ' . The novelty in this idea was 
that the traditional picture of Planck-length-sized additional spacelike dimensions 
(£p ~ 10 -33 cm) was abandoned, and the extra dimensions could have a size as large 
as 1 mm. The upper bound on the size of the proposed Large Extra Dimensions 
actually matched the smallest length scale down to which the force of the gravita- 
tional interactions, and thus their 1/r 2 dependence, had been measured. If extra 
dimensions of that size did exist, gravitational interactions would have a completely 
different dependence on r in scales smaller than 1 mm, however no gravitational 
experiment at that time could rule this out. 
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On the other hand, electromagnetic, weak and strong forces are indeed sensitive 
to the existence of extra dimensions. If, for example, gauge bosons were allowed 
to propagate in the extra-dimensional spacetime, their interactions would be mod- 
ified beyond any acceptable phenomenological limits unless the size of the extra 
dimensions was smaller than 10 _16 cm. This problem was resolved under the as- 
sumption that all particles experiencing this type of interactions, in other words, 
all ordinary matter, is restricted to live on a (3+l)-dimensional hypersurface, a 3- 
brane, that has a width along the extra dimensions of, at most, the above order. 
The 3-brane, playing the role of our four-dimensional world, is then embedded in 
the higher-dimensional spacetime, usually called the bulk, in which only gravity can 
propagate. 

In this model, the large volume of the extra dimensions can help us solve, or 
at least recast, the hierarchy problem: the traditional Planck scale, Mp, is only an 
effective energy scale derived from the fundamental higher-dimensional one, M*, 
through the following relation 1 

Mp~M* +n R n . (1) 

In the above, it has been assumed that each one of the extra spacelike, compact 
dimensions has the same size R. From the above, it becomes clear that, if the 
volume of the compact space, V <~ R n , is large, i.e if R » lp, then the (4 + n)- 
dimensional Planck mass, M», will be much lower than the 4-dimensional one, Mp. 
If one chooses M* = Mew, then the above expression provides a relation between 
the scale of gravity and the scale of particle interactions. 

The above idea makes use of geometrical features of the extra, compact space in 
order to connect two completely different energy scales. The same goal was achieved 
by an alternative, but similar idea, that was proposed a year later 3 (for some early 
works on brane-world models and their implications, see Refs. 4-10). According to 
this alternative proposal, the magnitude of the effective energy scale on the 3-branc, 
taken to be of the order of the Electro- Weak scale, follows from the fundamental 
higher-dimensional scale Mp after being suppressed by an exponential factor in- 
volving the distance of our observable brane from a hidden brane. In both models, 
however, the complete resolution of the hierarchy problem would demand also an 
explanation of why the volume of the compact space, or the inter-brane distance, 
has the value that leads to the observed ratio of Mp/Mew- 

In this review, we will concentrate on the scenario with Large Extra Dimen- 
sions. A number of experiments and theoretical analyses have tried over the years 
to put bounds on the size of extra, compact dimensions, and the produced bounds 
are constantly updated. In the regime KiJ, the extra dimensions 'open up' and 
Newton's law for the gravitational interactions is modified assuming a l/r™ +2 de- 
pendence on the radial separation between two massive particles. Torsion-balance 
experiments which measure the gravitational inverse-square law at short scales can 
provide limits on the size of the extra dimensions, or equivalently on the value of the 
fundamental scale M*. On the other hand, since gravitons can propagate both in 
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Tabic 1. Current limits on the size of extra, compact dimensions 



Type of Experiment / Analysis 




M, 


* > 




M, 


« > 




Collider limits on the production 

of real or virtual KK gravitons n. 12 . 13 


1.45 


TeV 


(n = 


2) 


0.6 TeV 


in 


= 6) 


Torsion-balance Experiments 14 


3.5 


TeV 


(n = 


2) 








Overclosure of the Universe 15 


8 


TeV 


(n = 


2) 








Supernovae cooling rate 16,17,18,19 


30 


TeV 


(n = 


2) 


2.5 TeV 


in 


= 3) 


Non-thermal production of KK modes 20 


35 


TeV 


(n = 


2) 


3 TeV 


in 


= 6) 


Diffuse gamma-ray background is. 21 . 22 


110 


TeV 


(n = 


2) 


5 TeV 


in 


= 3) 


Thermal production of KK modes 22 


167 


TeV 


(n = 


2) 


1.5 TeV 


in 


= 5) 


Neutron star core halo 23 


500 


TeV 


(n = 


2) 


30 TeV 


in 


= 3) 


Neutron star surface temperature 23 


1700 


TeV 


(n = 


2) 


60 TeV 


in 


= 3) 


BH absence in neutrino cosmic rays 24 










1-1.4 TeV 


in 


>5) 



the bulk and on the brane, massive Kaluza-Klein (KK) graviton states can modify 
both the cross sections of Standard Model particle interactions and astrophysical 
or cosmological processes. Assuming modifications that are below the current ob- 
servable limits puts bounds on the mass of the KK gravitons, and consequently, 
on the size of extra dimensions. In Table 1, we summarize some of those limits. 
The constraints from collider experiments, although more accurate, are particularly 
mild, while the cosmological and astrophysical ones are much more stringent; how- 
ever, they contain large systematic errors. If we ignore these errors, the latter type 
of constraints exclude by far even the n = 3 case, while low-gravity models with 
M* ~ 1 TeV are still allowed for n > 4. 

If indeed present, the extra dimensions will inevitably change our notion for the 
universe. The introduction of extra dimensions affects both gravitational interac- 
tions and particle physics phenomenology, and leads to modifications in standard 
cosmology. Already existing theories would need to be extended or modified in order 
to accommodate the effects resulting from the presence of extra dimensions (for an 
incomplete list of works on the cosmological and phcnomenological implications in 
theories with large extra dimensions, see Refs. 25-48). Similarly, the properties and 
physics of black holes are also bound to change in the context of a higher-dimensional 
theory. 

As in the four-dimensional case, it seems natural to assume that, when matter 
trapped on the brane undergoes gravitational collapse, a black hole is formed which 
is centered on the brane and extends along the extra dimensions. If the horizon of the 
formed black hole is much larger than the size of the extra dimensions, r# 3> R, the 
produced black hole is effectively a four-dimensional object. If, however, r# <C R, 
then this small black hole is virtually a higher-dimensional object that is completely 
submerged into the extra-dimensional spacetime. As we will see, these small black 
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holes have significantly modified properties compared to a four-dimensional black 
hole with exactly the same mass Mbh- for example, they are larger, colder and 
thus live longer compared to their four-dimensional analogues. Another striking 
consequence of the introduction of extra dimensions is that, by lowering the Planck 
scale M* closer to the Electro- Weak scale, the idea of the production of miniature 
black holes during high-energy scattering processes, with trans-Planckian center-of- 
mass energy s/s ^> M#, now becomes more realistic. Theoretical arguments have 
shown that the presence of the extra dimensions facilitates further the production 
of black holes during such collisions by increasing the production cross-section, 
thus leading to striking consequences for the high-energy interactions of elementary 
particles either at colliders or at cosmic rays. 

The produced black holes are characterized by a non-vanishing temperature 
Th, whose value is inversely proportional to the horizon radius. They decay by the 
emission of Hawking radiation, i.e. emission of elementary particles with rest mass 
smaller than T#. This is expected to be their most prominent observable signature 
with a characteristic thermal radiation spectrum and an almost blackbody profile. 
The non-trivial metric in the region exterior to the horizon of the black hole creates 
an effective potential barrier which backscatters a part of the outgoing radiation 
back into the black hole. The amount of radiation that finally reaches the observer at 
infinity depends on the energy of the emitted particle, its spin and the dimensionality 
of spacetime. The dependence on all the aforementioned parameters is encoded into 
the expression of a filtering function, the 'greybody factor' (lo), which is present 
in the radiation spectrum. The greybody factors can be important experimentally 
since they modify the spectrum in the low- and intermediate-energy regime, where 
most particles are produced, thus altering the characteristic spectrum by which we 
hope to identify a 'BH event'. In addition, as we will explain later in this review, 
by studying the Hawking radiation emitted by this type of small black holes, one 
would be able to 'read' the total number of dimensions that exist in nature. 

A higher-dimensional black hole emits radiation both in the bulk and on the 
brane. According to the assumptions of the theory with Large Extra Dimensions, 
only gravitons, and possibly scalar fields, can propagate in the bulk and thus, these 
are the only types of fields allowed to be emitted in the bulk during the Hawking 
evaporation phase. On the other hand, the emission on the brane can take the 
form of scalar Higgs particles, fermions and gauge bosons. From the perspective 
of the brane observer, the radiation emitted in the bulk will be a missing energy 
signal, while radiation on the brane may lead to experimental detection of Hawking 
radiation and thus of the production of small black holes. Nevertheless, in order to 
have a clear picture of the characteristics of the radiation spectrum on the brane, 
it is important to know exactly how much energy is lost in the bulk. 

As we mentioned above, the greybody factor depends on the dimensionality of 
spacetime; it also depends on whether the emitted particle is brane-localized or free 
to propagate in the bulk. The greybody factor is actually the outgoing transmis- 
sion cross-section associated with propagation in the aforementioned gravitational 
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background, however, due to the thermal character of the radiation spectrum, it is 
equal to the incoming absorption cross-section 49 . Therefore, all we need to do is to 
solve the equation of motion of a particle incident on the background metric that 
describes the black hole, either highcr-dimcnsional or four-dimensional. After the 
absorption coefficient is computed, the corresponding cross-section, and thus the 
greybody factor, can easily follow. 

We need to stress here that the above semiclassical calculation of Hawking emis- 
sion is only reliable when the energy of the emitted particle is small compared to 
the black hole mass, u) <C Mbh, since only in this case is it correct to neglect the 
back reaction to the metric during the emission process. This in turn requires that 
the Hawking temperature obeys the relation Th <C Mbh, which is equivalent to de- 
manding that the black hole mass Mbh S> M* . As the decay proceeds and the mass 
of the black hole decreases, inevitably this condition will break down during the fi- 
nal stages of the evaporation process. Nevertheless, for black holes of initial mass 
much larger than M* most of the evaporation process is within the semi-classical 
regime. 

We will start this article by reviewing, in Section 2, the existing literature on 
the creation of black holes during high-energy particle collisions, both in four- 
dimensional and higher-dimensional spacetimes. In Section 3, we will discuss the 
properties of the produced higher-dimensional black holes, namely the horizon ra- 
dius, the temperature and their life-time, and point out the differences between them 
and their four-dimensional analogues. The physics that governs the evaporation of 
these higher-dimensional objects, through the emission of Hawking radiation, is cov- 
ered in Section 4. Section 5 focuses on the emission of Hawking radiation directly on 
the brane: we start with the derivation of the master equation for the propagation 
of fields with arbitrary spin in the induced-on-the-brane black hole background, 
and then we present all existing results in the literature for the emission of scalars, 
fcrmions and gauge bosons during the spin-down and Schwarzschild phases of the 
life of the black hole; both analytical and numerical results for the greybody factors 
and radiation spectra are presented as well as exact results for the number and 
type of fields emitted on the brane as a function of the dimensionality of space- 
time. Section 6 deals with the emission of Hawking radiation in the bulk and the 
question of the amount of the missing energy: analytical and numerical results on 
the emission of bulk scalar fields, including greybody factors and radiation spectra, 
are reviewed, and the ratio of the missing energy over the 'visible' one emitted on 
the brane is presented for the Schwarzschild phase and for different values of the 
number of extra dimensions. Our conclusions are summarized in Section 7. 

2. High-energy Collisions and Black Hole Creation 

If extra dimensions do exist and the fundamental scale of gravity is much lower 
than the traditional Planck scale Mp, then in the near future we will witness col- 
lisions of particles with trans-Planckian energies, i.e. energies larger than the fun- 
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damental Planck scale M*. Although, as we approach M», quantum gravity effects 
become important for elementary particles, another semi-classical regime opens up 
for center-of-mass energies of the colliding particles \fs >• M* ■ the products of such 
a collision would have a mass much larger than the scale of quantum gravity and 
for these objects (as for any macroscopic object) quantum gravity effects can be 
safely ignored. 

The high-energy scattering of particles and the nature of the products of such 
collisions have been investigated in the framework of the General Theory of Relativ- 
ity 50,51,52,53,54^ string Theory 55 < 56 and Quantum Gravity 57 ' 58 . The usual approach 
followed is that the relativistic colliding particles (or black holes) may be described 
by two Aichclburg-Sexl gravitational shock waves. In the limit of moving velocity 
equal to the speed of light, the two particles can be considered massless and the cur- 
vature is zero except on the null plane of their trajectory. If the impact parameter b 
is larger than the Schwarzschild radius rn, that corresponds to the center-of-mass 
energy of the two particles, elastic and inelastic processes will in general take place 
between the two shock-waves accompanied by the exchange of gravitons. If on the 
other hand, b < m, then, according to the General Theory of Relativity and the 
"hoop-conjecture" 59 , strong gravitational effects will dominate and a black hole 
will be formed. The production cross-section, in this high-energy limit, is given by 
the geometrical cross-section 



i.e. by the "target" area defined by the impact parameter. 

As reported in Ref. 54, Penrose found a lower bound on the mass of the black hole 
produced during the collision of two particles moving at the speed of light 60 . On 
the union of the two null planes, that describe the trajectories of the two particles, 
an apparent horizon is formed, with an area of 327r/i 2 where /i is the energy of each 
particle in the center-of-mass coordinate frame. This sets a lower bound on the area 
of the event horizon (that according to the Cosmic Censorship hypothesis 61 lies 
outside the apparent horizon) and thus on the mass of the produced black hole; this 
is found to be 



leading to the conclusion that at least 71% of the initial center-of-mass energy of the 
colliding particles will be bound into the black hole. The perturbative analysis done 
subsequently in Ref. 54 determined the amount of energy spent into gravitational 
radiation during the collision to be 16% of the total energy, thus raising the estimate 
of the mass of the black hole to 84% of y/s. The above results were derived under 
the assumption of head-on collision (b = 0) in four-dimensional spacetime. An 
alternative study of the high-energy collision of a particle with a Schwarzschild 
black hole in four dimensions 62 led to a similar result: the emission of gravitational 
radiation is approximately 13% of the initial energy of the system; in the case of 



a ~ nb 2 ~ nrj^ , 



(2) 



M BH > 4= 2 M 



(3) 
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a rotating black hole, however, this percentage can become as large as 35%, but 
decreases as the impact parameter increases. 

After the theories with large extra dimensions and a low-scale gravity were 
proposed 1 , the idea of the creation of black holes from the collision of particles with 
trans-Planckian energies was revived 63 . The need for the update of the above lower 
limit on the mass of the black hole for an arbitrary number of spacelike dimensions 
was obvious. Moreover, realistic calculations of the rate of production of black holes 
during particle collisions demanded the generalization of the same analysis for non- 
zero impact parameters. By investigating the formation of a closed-trapped surface 
around the colliding particles/shock- waves, the authors of Ref. 64 have shown that 
such surfaces do indeed form for both D = 4 and D > 4 upon particular assumptions 
for the value of the impact parameter. Of particular importance was their result 
that, in D = 4, black holes form if, and only if, b < b max ~ 0.8rn, a result which 
reduces the value of the cross-section to a ~ 0.65"7rr^, and sets the range of the 
black hole mass to 

Mbh — (0.71 — 0.45)\/s , for b={0,b max }. (4) 

Their higher-dimensional analysis, performed only for head-on collisions (b = 0), 
revealed that the lower bound on the mass of the produced black hole decreases 
also with the number of extra dimensions reaching 0.58-y/s for D = 11. In Ref. 65, 
it was found that, for head-on collisions and as D increases, the circumference of 
the region into which the mass M must be compacted to produce a black hole is 
a decreasing fraction of 2nrn- The same authors also showed that, for 6^0, the 
value of b max that can produce a black hole increases with D and is given by: 

b max ~2-^ D -Vr H . (5) 

This leads to the enhancement of the production cross-section, nevertheless, the 
mass of the produced black hole will still be suppressed for large values of D, in 
agreement with the results of Ref. 64. A recent numerical analysis 66 has produced 
some new estimates for the energy lost in the form of gravitational radiation during 
a head-on collision in a higher-dimensional spacetime: 13% for D = 4 to 8% for 
D = 10. While the former result is in agreement with Refs. 54 and 62, the latter 
seems to disagree with the results produced by the authors of Refs. 64 and 65. The 
apparent disagreement can be resolved only under the assumption that a significant 
part of the energy lost in the collision has a form different from that of gravitational 
radiation. 

The formation of closed-trapped surfaces was studied analytically in Ref. 67 for 
the more realistic case of finite-front shock-waves, and was shown to indeed take 
place for D > 4 and for an arbitrary impact parameter. The effect of the angular 
momentum was looked at in Ref. 68, where it was argued that the production cross- 
section gets further enhanced when the spinning of the black hole is taken into 
account. Finally, the evolution of the closed-trapped surfaces in time was studied 
in Ref. 69. 
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But where and when these black holes may be produced? In the context of theo- 
ries with large extra dimensions, the trans-Planckian energy regime may lie slightly 
above the TeV scale. This raises the exciting possibility that particle collisions with 
trans-Planckian energies may take place, in the near future, at ground-based accel- 
erators 70 > 71 > 72 > 73 > 74 j or they may even already take place in the atmosphere of the 
earth 75 > 76 . 77 . 78 . T n the first case, accelerated particles (protons or nuclei) collide 
at center-of-mass energies \fs > M* at ground-based accelerators; if the impact 
parameter is smaller than a critical value b c , a higher-dimensional black hole will 
be formed. The same holds in the case where highly-energetic cosmic ray particles 
with energies up to 100 TeV collide with particles in the atmosphere of the earth. 
Neutrinos being scattered by nuclei are expected to be the most effective source of 
black hole production from cosmic ray particles due to the absence of any QCD-type 
contaminating effects and small Standard-Model cross-sections. 

The corresponding black hole production cross-sections are calculated in both 
cases in the same way: first, any pair of partons that pass within the 

Schwarzschild radius, i.e. have b < b c ~ ^(s), can lead to the production of a 
black hole with cross-section ct^ ~ wr%(s), where y/s is the center-of-mass energy 
of the colliding particles. However, in a realistic collision, the colliding particles con- 
sist of more than one partons, therefore, a summation must be made over all parton 
pairs that carry enough energy to produce a black hole of a minimum mass M min . 
For this reason, the fraction of the center-of-mass energy that each parton carries 
must be taken into account via the use of the parton distribution functions fi(x), 
that give the probability of finding a parton with a fraction x of the momentum 
of the colliding particle. The total cross-section for the production of a black hole 
from the collision of two particles A and B may then be written as 70 ' 71 



where r — XiXj is the parton-parton center-of-mass energy squared fraction, and 
y/T m s is the minimum center-of-mass energy necessary for the creation of the min- 
imum black hole mass, M min ~ y/r m s. In the above, it has been assumed that 
both colliding particles are composite - in case one of the particles is elementary 
(i.e. neutrinos) this expression is simplified 76 . Early estimates found particularly 
large production cross-sections of black holes: high-energy neutrinos scattered by 
nuclei at the earth's atmosphere give a production cross-section, for a black hole 
with M min = M* = 1 TeV, which is two orders of magnitude larger than the cross- 
section of any similar Standard Model process 76 . On the other hand, proton- proton 
scattering at the LHC with M. t — 1 TeV was found to lead to a production cross- 
section of 10 5 fb for Mmin = 5 TeV, a large cross-section by new physics standards; 
for M min — 10 TeV, the cross-section reduces to 10 fb. These estimates gave a 
substantial boost towards the further study of the production and phenomenology 
of black holes in theories with extra dimensions (see Refs. 79-106). 

We should note here that the above production cross-sections decrease if one 
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takes into account the results of the studies on the formation of closed-trapped 
surfaces mentioned earlier in this section. By taking b c = m and not b c = b max , 
where b max is only a fraction of the Schwarzschild horizon, we obviously obtain an 
overestimate of the production cross-section. Moreover, assuming that the available 
energy for the production of black hole is the whole of the center-of-mass energy y/s 
is again an over-simplification since part of that energy is bound to be lost in the 
form of gravitational radiation during the collision 54 . A claim for an additional sup- 
pression of this process was made 107 according to which an exponential suppression 
factor involving the Euclidean action of the system should be included. However fur- 
ther studies 72 ' 67 ' 108 argued that the black hole creation process from the collision 
of two particles was not classically forbidden and this suppression factor should 
not be taken into account. Subsequent analyses 109 > 110 calculated quantum correc- 
tions to the semi-classical cross-section and showed that these are indeed small. 
From a different perspective, the authors of Ref. Ill have claimed that the use of 
the generalized uncertainty principle leads to a radical increase in the minimum 
amount of energy that is necessary for the creation of a black hole, thus render- 
ing unlikely the production of black holes at next-generation colliders such as the 
Large Hadron Collider (LHC) with a center-of-mass energy of 14 TeV. Recently, 
the validity of the description of the creation of a black hole from the collision of a 
pair of Aichelburg-Sexl shock-waves was also questioned 112 , and the argument that 
strong-curvature and quantum gravity effects may significantly alter the geometric 
cross-section estimate was put forward. 



3. Properties of the Higher-Dimensional Mini Black Holes 

As we will see below, in the absence of extra dimensions, the creation of a semi- 
classical black hole demands extremely large center-of-mass energies, which are far 
beyond our technical abilities. Reducing the mass of the black hole down to accessi- 
ble energy scales simply leads to an unnaturally small Schwarzschild radius, which 
can never be attained. On the other hand, the presence of the extra dimensions 
facilitate the creation of black holes since it lowers the scale of quantum gravity, 
thus allowing the production of semi-classical black holes at lower energies, and 
increases the corresponding Schwarzschild radius for a given center-of-mass energy, 
thus making the black hole creation regime, b < ru, more easily accessible. 

Of particular importance, and simplicity, are the higher-dimensional black holes 
that have a horizon radius much smaller than the size of the extra dimensions, 
tr *C R. In this case, these mini black holes are completely submerged into a D- 
dimensional spacetime that, to a very good approximation, has one timelike and 
D — 1 non-compact spacelike coordinates. If we further assume that the produced 
black hole is spherically-symmetric, i.e. non-rotating, the gravitational background 
around this black hole is given by a generalized Schwarzschild line-element of the 
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form 113 

V r ) 



ds z = 



1 _(2L) n+1 ] ' 1 ; 



dr 2 - r 2 dQ 2 2+n , (7) 



where n stands for the number of extra, spacelike dimensions that exist in nature 
(D = 4 + n) , and dQ\ +n is the area of the (2 + n)-dimensional unit sphere given by 

dQ 2 2+n = d9 2 n+1 + sin 2 9 n+1 (d6l + sm 2 6 n (... + sm 2 6 2 (c$ 2 + sin 2 <9i dip 2 ) ...U .(8) 



In the above, < tp < 2ir and < 6i < ir, for i = 1, n + 1. The line-element (7) 
can be easily shown to satisfy the vacuum (4 + n)-dimcnsional Einstein's equations. 
The black hole is assumed to be bound on a 3-brane, our four-dimensional world, 
nevertheless, the tension of the brane is assumed to be much smaller than the black 
hole mass and thus it can be ignored in our analysis. 

By using an analogous approach to the usual 4-dimensional Schwarzschild calcu- 
lation, i.e. by applying Gauss' law in the (4 + n)-dimensional spacetime, we obtain 
the following relation between the horizon radius and the mass Mbh of the black 
hole 113 

rfl = J_«*[Hffl)f (9 ) 
H ^M, \ M, ) \ n + 2 J [ ' 

We notice that, for ji^O, the relation between rjj and Mbh is not linear anymore, 
and that it is the fundamental Planck scale M* that appears in the expression of 
the horizon radius and not the four-dimensional one Mp. The latter feature is the 
main reason for the fact that extra dimensions facilitate the creation of low-mass 
black holes, as we will shortly see. 

Before elaborating on this last point, we need to make another comment first: 
in order to be able to ignore quantum corrections in our calculations and study 
the produced black holes by using semi-classical methods, the mass of the black 
hole must be, at least, a few times larger than the scale of quantum gravity M*. 
Therefore, if we assume that M* = 1 TeV, a safe limit for the mass of the produced 
black hole would be 70 ' 73 Mbh — 5 TcV. By keeping fixed the mass of the produced 
black hole, we may calculate the value of the horizon radius as a function of n; these 
values are given in Table 2. 

Table 2. Black hole horizon radii for different values of n 



n 1 2 


3 


4 


5 


6 


7 


r H (1(T 4 fm) 4.06 2.63 


2.22 


2.07 


2.00 


1.99 


1.99 



One may easily conclude, from the above, that during the collision of two par- 
ticles, with a center-of-mass energy y/s > 5 TeV, a black hole may be formed if 
the particles pass within an area of radius 10~ 4 fm; this is merely a sub-nuclear 
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distance attainable at particle physics experiments. On the other hand, in the ab- 
sence of extra dimensions, the lightest semi-classical black hole would have a mass 
of, at least, a few times the four-dimensional Planck mass, Mp ~ 10 16 TeV; the 
creation of a black hole from a high-energy collision would then demand center- 
of-mass energies higher than Mp, a requirement which is far beyond the reach of 
any present and future accelerator. Overlooking for the moment the fact that a 
four-dimensional black hole with Mbh < Mp would not be a classical object, we 
may ask what would be the value of the Schwarzschild radius for an object with 
mass M = 5 TeV being produced in four dimensions in a high-energy collision. We 
find that the Schwarzschild radius for such an object has the extraordinary value 
of tr = 1-3 x 1CP 50 m, i.e. 35 orders of magnitude smaller than the radius of the 
proton. 

The modified properties of a higher-dimensional, Schwarzschild-like black hole, 
compared to those of a four-dimensional one with the same mass, were first studied 
in Ref. 114. The fact that the Schwarzschild radius in D > 4 dimensions is larger 
than the one in D = 4, for a given mass Mbh, was first pointed out in there, and 
further implications of the existence of extra dimensions on the temperature, life- 
time and entropy of the black hole were investigated. We will now stop and look 
in some detail at the temperature of a higher-dimensional black hole to see where 
the difference from the four-dimensional case lies. The temperature of a (4 + n)- 
dimensional black hole is given by the expression 113 

T H = ^. (10) 

47T TH 

Let us assume again that M* = 1 TeV, and that the produced black hole has a 
mass Mbh = 5 TeV. By using Eq. (10) and the entries of Table 2, we may easily 
calculate the temperature of the produced black hole for different values of n; the 
results arc displayed in Table 3. 



Tabic 3. Black hole temperatures for different values of n 



n 


1 


2 


3 


4 


5 


6 7 


T H (GeV) 




179 


282 


379 


470 


553 629 



As we will shortly see, a black hole with a temperature Th emits thermal radia- 
tion, the so-called Hawking radiation 115 , through the emission of ordinary particles. 
This leads to the decay of the black hole and finally to its evaporation. The radia- 
tion spectrum has a blackbody profile with the peak of the curve being at energies 
very close to its temperature. The temperature values displayed in Table 3 all lie 
in the GeV regime, which is the energy range that present and next-generation 
experiments can probe. We may, thus, conclude that the presence of extra dimen- 
sions, not only facilitates the creation of a black hole at a high-energy collision, but 
also renders more likely the detection of their most prominent feature, the emitted 
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Hawking radiation. 

We should note here that Hawking radiation emitted by larger, and thus effec- 
tively four-dimensional, astrophysical black holes has never been observed. Taking 
Mbh — 3M Q , which is the current upper limit on the mass of neutron stars, and 
using Eq. (10) with n — 0, we find that Th — 10 -12 eV = 20 nK; this is a very low 
temperature corresponding to a very low energy frequency, or a very large wave- 
length, that unfortunately cannot be detected. In the absence of extra dimensions, 
the only black holes that could emit radiation today at detectable frequencies, are 
primordial ones that were created at the early universe and can have a much lower 
mass. For example, a black hole with a mass Mbh = 10 15 gr should give a Hawking 
radiation spectrum with a peak in the area of 10-100 MeV; still, such a radiation 
has not been yet observed. 

Combining Eqs. (9) and (10), we may conclude that, since a (4 + n)-dimensional 
small black hole has a horizon radius much larger than a four-dimensional one with 
the same mass, it will have a temperature in D > 4 which is much lower than the 
one in D = 4 114 . Indeed, if a 'quantum' black hole with Mbh = 5 TeV had been 
allowed to exist in nature in D = 4, its temperature would have been 30 orders 
of magnitude larger than the entries in Table 3. The larger the temperature of 
the black hole is, the faster its decay rate is - through the emission of Hawking 
radiation - and thus the shorter its lifetime. Astrophysical black holes with masses 
Mbh > 3M emit radiation with an extremely small rate a and their lifetime, given 
by the four-dimensional relation 



turns out to be much larger than the age of the universe. On the other hand, the 
same formula tells us that tiny black holes (as the ones that might have been created 
in the early universe in D = 4) have an extremely short life-time, due to their huge 
emission rate, and must have decayed and evaporated long time ago. However, for 
small black holes, the presence of extra dimensions modifies also their lifetime, and 
it is now given by 113 



The appearance of the low energy scale M* in Eq. (12), instead of Mp as in Eq. (11), 
leads to black hole lifetimes much longer than the one in four dimensions. Neverthe- 
less, for black holes with masses in the area of a few TeV, the lifetime is still a tiny 
fraction of a second: for Mbh — 5 TeV, the lifetime ranges from 1.7 x 10~ 26 sec (for 
n = 1) to 0.5 x 10~ 26 sec (for n = 7), while for Mbh = 10 TeV, the corresponding 
lifetime interval is from 1.6 x 10~ 26 sec (for n = 1) to 1.2 x 10~ 26 sec (for n = 7). 




(11) 




(12) 



a In addition, those black holes have a temperature much smaller than the one of the Cosmic 
Microwave Background Radiation, T = 2.73 K = 2.3 X 10~ 4 eV, therefore, they actually absorb 
radiation from their environment instead of emitting. 
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Despite their short lifetime, the fact that these small black holes might be created 
during collisions on, or close to, the Earth's surface makes the possibility of their 
observation a more realistic prospect. 



4. Evaporation of Higher-Dimensional Black Holes 

The emission of Hawking radiation is in fact compatible with the well-known result 
of the General Theory of Relativity that nothing can escape from inside the horizon 
of a black hole. The Hawking radiation can be conceived as the creation of a virtual 
pair of particles just outside the horizon of the black hole: the particle with the 
positive energy escapes to infinity while the antiparticle (the one with the negative 
energy) falls into the BH, where it can exist as an ordinary particle. The spectrum of 
the Hawking radiation coming from a black hole with temperature Th is a thermal 
one with an almost blackbody profile. The flux spectrum, i.e. the number of particles 
emitted per unit time, from a higher-dimensional spherically-symmetric black hole of 
the type (7), can be easily found by generalizing the corresponding four-dimensional 
expression 115 for a higher number of dimensions. It is given by 

dt aj ' n{ ' exp (u/T H ) ± 1 (2tt)™+ 3 ' { ' 

where s is the spin of the emitted degree of freedom and j its angular momentum 
quantum number. The spin statistics factor in the denominator is —1 for bosons and 
+1 for fermions. For massless particles, \k\ = uj and the phase-space integral reduces 
to an integral over the energy of the emitted particle u>. For massive particles, 
\k\ 2 = uj 2 — to 2 , and the energy in the denominator now includes the rest mass of 
the particle: this means that a black hole temperature Th > to is necessary for 
the emission of a particle with mass to. We should note here that, as the decay 
progresses, the black hole mass decreases and the Hawking temperature rises. It is 
usually assumed that a quasi-stationary approach to the decay is valid - that is, 
the black hole has time to come into equilibrium at each new temperature before 
the next particle is emitted. We will make this assumption also here. 

The power spectrum, i.e. the energy emitted per unit time by the black hole, 
can be easily found by combining the number of particles emitted with the amount 
of energy they carry. It is given by 

dE^juj) v (s) uj rf"+ 3 fc 

dt a j,n\ UJ ) exp (u/Th) ± 1 (27r) n + 3 ' 1 ' 

Both expressions, Eqs. (13) and (14), contain an additional factor, <t^(w), which 
does not usually exist in a typical blackbody spectrum, or, more accurately, is just 
a constant standing for the area of the emitting body. In contrast, this factor here 
depends on the energy of the emitted particle, its spin and its angular momentum 
number. It may therefore significantly modify the spectrum of the emitted radiation 
and, for that reason, is called the 'greybody' factor. Equally important is the fact 
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that this coefficient depends also on the number of extra dimensions and therefore 
encodes valuable information for the structure of the spacetime around the black 
hole including the dimensionality of spacetime. 

The distortion of the blackbody spectrum, or in other words, the presence of 
the greybody factor in the radiation spectrum of a black hole, can be attributed 
to the following fact: any particle emitted by a black hole has to traverse a strong 
gravitational background before reaching the observer at infinity, unlike to what 
happens with a black body in flat spacetime. The radiation spectrum is bound to 
depend on the energy of the propagating particle and the shape of the gravitational 
barrier as these are the parameters that will determine the number of particles that 
manage to reach infinity. In order to illustrate the above, let us assume, for example, 
that a scalar field of the form <fi (t, r, 9i,(p) — e~ lut R w t(r) Ye(il), where Ye(SY) is the 
(4 + n)-dimensional generalization of the usual spherical harmonic functions 116 , 
propagates in the background of Eq. (7). Its radial equation of motion may then be 
written in a Schr6dinger-like form 



A 2 
dy 2 



-c 2 + 



l)h(r) 



R u e(y) = 0, (15) 



in terms of the 'tortoise' coordinate 

The quantity inside square brackets in Eq. (15) gives the effective potential barrier in 
the area outside the horizon of the black hole. It clearly depends on all parameters 
mentioned above: u), j (being equal to the orbital angular momentum number I 
in the case of a scalar particle) and n. As a similar analysis will shortly show, it 
also depends on the value of the spin s, for a non-zero spin particle. By simple 
inspection, one may see that the barrier lowers for larger energy w, while it rises 
for higher angular momenta I - the dependence, however, on the number of extra 
dimensions is more subtle. 

The greybody factor a^^uj) then stands for the corresponding transmission 
cross-section for a particle propagating in the aforementioned background. This 
quantity can be determined by solving the equation of motion of a given particle and 
computing the corresponding absorption coefficient Aj . Then, we may write 117 

(s) 2» 7 r("+ 1 >/ 2 r[(n+l)/2] (2j + n + 1) (j + n)! (s) 2 

°j,nW- Jl \ A j I ' y ll > 

It will be useful to rewrite the above expression for the greybody factor as 

where Nj is the multiplicity of states corresponding to the same partial wave j, 
given for a (4 + n)-dimensional spacetime by 118 
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and Ah is the horizon area of the (4 + n)-dimensional black hole defined as 

A H = r H +2 / d V TT / sin fe k+1 d9 k+1 
Jo fe=1 Jo 



n+l 



- r* (27T) JJ_ V- r[(fc + 2)/2] 



^+ 2 (27r)7r(" +1 )/ 2 rf^) \ (20) 



From Eq. (18), we may see that the greybody factor is indeed proportional to 
the area of the emitting body, as in the case of blackbody emission, nevertheless 
additional factors change, in principle, this simple relation by adding an explicit 
dependence on u, rn, j and n. We should also note that a^{u) has the same 
dimensions as Ah, therefore its dimensionality changes as the number of extra 
dimensions n varies. Normalizing its expression to the area of the horizon of the 
(4 + ra)-dimensional black hole would therefore be useful in comparing its values for 
different n. 

Equation (14), for energy emission in the higher-dimensional spacetime, may 
take a simpler form when written in terms of the absorption coefficient. Then, it 
reads 

dE^jw) ^ ()2 w du 

dt 2^ ^ I A J I exp (cj/T B h) ± 1 2tt ' { ' 

In the above, we have assumed that the emitted particles are massless, an assump- 
tion that will be made throughout this review. This simple, alternative form of 
the power spectrum has led part of the community to referring to the absorption 
probability |^4j S ' ) | 2 as the greybody factor, this however will not be followed here. 

As wc mentioned above the greybody factor modifies the spectrum of emitted 
particles from that of a perfect thermal blackbody 115 . Four-dimensional analyses 
for Schwarzschild 119 . 12 o,i2i anc j Kerr 122 black holes have determined, both an- 
alytically and numerically, the greybody factors for particles of different spin. In 
the simplest case of a non-rotating, Schwarzschild black hole, geometric arguments 
show that, in the limit of high energy uj (the geometrical optics limit), Xj crj^(w) 
is a constant independent of ui 123 (for more details, see Sections 5 and 6). In that 
case, the spectrum is exactly like that of a blackbody for every particle species inde- 
pendently of their spin s. The low-energy behaviour, on the other hand, is strongly 
spin-dependent and energy-dependent, and the greybody factors are significantly 
different from the geometrical optics value 119 > 12 M24 ^he result is that both the 
power and flux spectra peak at higher energies than those for a blackbody with the 
same temperature. Finally, the spin dependence of the greybody factors means that 
they are necessary to determine the relative emissivities of different particle types 
from a black hole. 
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For a charged black hole, early works 125 > 126 > 127 showed that the black hole will 
quickly discharge through a Schwinger-type pair-production process. Moreover, the 
charge of the black hole could affect the geometry of spacetime and thus the emission 
of uncharged particles only for black holes with masses larger than 10 5 Mq. A few 
years later, Page 128 showed that during the evaporation phase, the charge of the 
black hole fluctuates. The electrostatic potential between particle and antiparticle, 
together with the charge fluctuations, reduces the emitted flux and power of charged 
particles, bringing it down to 7% of the emission rate for similar but uncharged 
particles. In the same work, the effect of the rest mass m of the emitted particle 
was studied and found that it also reduces the emission rate: for electrons and 
muons, the reduction could be as large as 50% for black hole masses of 0(1O 16 gr). 

In the light of the recent speculations of the possible creation of higher- 
dimensional black holes during high-energy particle collisions, the generalization 
of the existing four-dimensional analyses for the calculation of greybody factors to 
a higher number of dimensions becomes imperative. On the one hand, the existence 
of extra dimensions might significantly change all that we know about the emis- 
sion of particles from a black hole, including the rate of emission and the type of 
particles emitted. Moreover, if the dependence of the radiation spectrum on the 
dimensionality of spacetime is strong enough, then, a possible detection of such a 
spectrum can give us valuable information on the dimensionality of spacetime. 

If, eventually, small higher-dimensional black holes may be created in high- 
energy particle collisions, then according to Rcfs. 70 and 73, the produced black 
holes will go through a number of phases before completely evaporating. These are: 

• The balding phase : the black hole emits mainly gravitational radiation and sheds 
the 'hair' inherited from the original particles, and the asymmetry due to the 
violent production process. 

• The spin-down phase : the typically non-zero impact parameter of the colliding 
partons leads to black holes with some angular momentum about an axis perpen- 
dicular to the plane. During this phase, the black hole loses its angular momentum 
through the emission of Hawking radiation and, possibly, through superradiance. 

• The Schwarzschild phase : a spherically-symmetric black hole loses energy due to 
the emission of Hawking radiation. This results in the gradual decrease of its mass 
and the increase of its temperature. 

• The Planck phase: the mass and/or the Hawking temperature approach - a 
theory of quantum gravity is necessary to study this phase in detail. 

An important question related to the aforementioned life stages of the black 
hole is how much energy is spent during each one of those phases. We have already 
seen how the initial four-dimensional estimate of 16% of the total energy emitted 
in the form of gravitational radiation during the balding phase, can go up to 55% 
depending on the value of the impact parameter and the dimensionality of space- 
time. Crude estimates of the corresponding percentages in four dimensions give 25% 
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and 60% for the spin-down and Schwarzschild phase 70 > 73 j respectively. However, it 
becomes clear that exact analyses performed in the higher-dimensional spacetime, 
in conjunction with the results derived already for the balding phase, can radically 
change those numbers. 

In what follows, we will concentrate on the second and third phase of the life of 
the black hole, that is the spin-down and Schwarzschild phase. We will ignore the 
charge of the black hole: although we expect the aforementioned 7% figure, for the 
emission rate of charged particles, to change with the addition of extra dimensions, 
we anticipate that any increase would still allow us to ignore the emission of charged 
particles, at least at first approximation. We will also study the emission of only 
massless particles: the mass of the black hole, being in the trans-Planckian regime 
of a few TeV's, leads to temperatures that are much larger than the rest masses of 
all known particles (see Table 3). 

Another comment is in order here: the presence of the brane has so far been 
ignored under the assumption that its tension is much smaller than the mass of the 
black hole and thus it does not change the gravitational background. Nevertheless, 
the brane plays another role, that of our four-dimensional world where, according to 
the assumptions of the Theories with Large Extra Dimensions, all ordinary particles 
(fermions, gauge bosons and Higgs fields) are localized. A (4 + n)-dimensional black 
hole emits Hawking radiation both in the bulk and on the brane. Since only gravi- 
tons, and possibly scalar fields, live in the bulk, these are the only particles that can 
be emitted in the bulk. On the other hand, the emission of brane-localized modes 
include zero-mode scalars, fermions, gauge bosons and zero-mode gravitons. From 
the observational point of view, it is much more interesting to study the emission of 
brane-localized modes, however, the emission in the bulk is equally important since 
it provides answers to questions like how much energy is available for emission of 
Standard Model particles on the brane. 

The brane-localized modes and bulk modes live in spacetimes of different dimen- 
sionality. Whereas the bulk modes have access to the whole of the (4+n)-dimcnsional 
spacetime, the brane modes live in a four-dimensional slice of it, which is the pro- 
jection of the higher-dimensional spacetime on the brane. In order to study the 
emission of both types of modes, we need to write down the corresponding equa- 
tions of motion of those fields. In the next section, we will first derive the master 
equation describing the motion of a field with arbitrary spin s in the background 
of a higher-dimensional Kerr-like black hole projected onto our brane; then, we will 
move on to review the so far derived results for the emission of Hawking radiation 
from rotating and non-rotating, uncharged (4 + n)-dimensional black holes on the 
brane. A similar review of the emission of Hawking radiation in the bulk is left for 
Section 6. 
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5. Master Equation and Hawking radiation on the Brane 

The gravitational background around a (4 + n)-dimensional, rotating, uncharged 
black hole was found by Myers & Perry, and is given by the following line-element 113 

ds 2 = (l ^V* 2 + t^^dtdv- ^dr 2 -ZdO 2 

_ ( r 2 +a 2 + °^ Sin ' 6 ) sin 2 6 d^ 2 - r 2 cos 2 6 dQ n , (22) 



£) r n-l 



where 



A = r 2 + a 2 --^ T , E = r 2 +a 2 cos 2 , (23) 



and <if2„ is the line-element on a unit n-sphere. The mass and angular momentum 
(transverse to the r^-plane) of the black hole are given by 

with G being the (4 + n)-dimensional Newton's constant, and ^4„+2 the area of a 
(n + 2)-dimensional unit sphere given by 

2n (n+3)/2 

An+2 = T[(n + 3)/2]- (25) 

We will assume that the line-element (22) describes successfully a small, higher- 
dimensional black hole during its spin-down phase. At the end of this phase, the 
angular momentum will be lost, and the same line-element, in the limit a — > 0, 
will describe the following spherically-symmetric Schwarzschild phase of the life of 
the black hole. As mentioned earlier, in this section, we will be interested in the 
emission of Hawking radiation on the brane, therefore we need to determine the 
line-element in which the brane-localized modes propagate. This can be found by 
fixing the values of the additional angular coordinates describing the compact n 
dimensions, and, then, the induced-on-the-brane line-element is described again by 
Eq. (22) but with the last term omitted. 

In order to derive a master equation describing the motion of a field with arbi- 
trary spin s in the aforementioned projected background, we need to make use of 
the Newman-Penrose formalism 129 - 130 that puts spinor calculus and gravitational 
quantities in the same framework. We first need to choose a tetrad basis of null 
vectors n^, to m , m**), where I and n are real vectors, and m and rh are a pair of 
complex conjugate vectors. They satisfy the relations: 1 • n = 1, m • m = —1, with 
all other products being zero. Such a tetrad basis is given below: 

^ = T 2+a2 1 -) n"= r 2 + ° 2 -— —) 

m/ 1 = f iasinO, 0, 1, -^—r) —=- , mt 1 = (-iasm6, 0, 1, -^-) , (26) 

V sm0/ V2p v sin0/ V2p* 
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where p = r + iacosO. The X a b c coefficients, which are used to construct the spin 
coefficients, are defined as: 



A q6c = {e b )i,j {e a ) l (e c ) J - {e a ) J (e c y 



(27) 



where e a stands for each one of the null vectors and (i,j) denote the components of 
each vector. To determine X a bc we will also need the expressions of the null vectors 
with their index down, and these are: 



(l, ~, 0, -asm 2 9) 



f A 1 a sin 2 9 A 



m,. 



iasin9, 0, — S, — isin9(r +a) 



V2S ' 2 
1 



0, 



2S 



V2(> 



(28) 



Then, the non-vanishing components of the X a bc coefficients are found to be: 



A 122 = -^(|c? r A-rA), A 
\/2a 2 cos 9 sin 9 



132 



sin 9 



Al43 — - , 
P 



A213 — — 
A324 = — 



ia cos 9 A 

m 2 ' 



A 2iacos9 
A243 = _ 2^' Asi4 = - 



A334 — 



r cos 9 + ia 
V2p 2 sin 



(29) 



The above components must be supplemented by those that follow from the sym- 
metry X a b c = — A c ba and the complex conjugates obtained by replacing an index 3 
by 4 (or vice versa) or interchanging 3 and 4 (when they are both present). 

We may now compute the spin coefficients defined by j a b c = (X a bc + A ca h — 
A&ca)/2. Particular components, or combinations, of the spin coefficients can be 
directly used in the field equations 129 - 130 . These are found to have the values: 



k = a = A = v 
iasm9 



7T 



7 : 



0, 

&.A 



A 



2Y,p* ' 
cot 6> 



a = it - [3* . 



V2> 2 ' 4S 2£p* ' " 2V2p 

In what follows, we will also employ the Newman-Penrose operators: 

b 



(30) 



(r 2 +a 2 ) 8 d_ a_d_ 
A dt + dr + A dip ' 



(r 2 + a 2 ) 9 A d o_ d_ 
2S di ~ 2E + 2E dip ' 



(5 = 



1 



ia smtf — + 



d 



y^2p V dt 86 sin0 9^ ' 

and make use of the following field factorization: 

tf s (t, r, 6, V ) = e-*"* e mip R s (r) (0) , 



(31) 
(32) 
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where S™J8) are the so-called spin- weighted spheroidal harmonics. For s = 0, these 



angular cigcnfunctions reduce to the spheroidal harmonics 131 , while for au; = 0, 
they take the form of the spin- weighted spherical harmonics 132 . In what follows, we 
will address separately the equation of motion of each type of field (gauge bosons, 
fermions and scalars) in the background of the projected, Kerr-like black hole. 

Gauge Bosons (,s = 1). In the Newman-Penrose formalism, there are only three 
'degrees of freedom' for a gauge field, namely $o = -F13, $1 = (F12 + F^)/2 and 
$2 = -F42, in terms of which the different components of the Yang-Mills equation 
for a massless gauge field are written as: 



(D - 2p) $1 - (5* + 7T - 2a) $o=0, 
(5 - 2t) $! - (A + n - 2 7 ) $o = 0, 
(D-p) $2 - (5* + 2ir) $x =0, 
(S-t + 2(3) $2 - (A + 2/x) *i=0, 



(33) 
(34) 
(35) 
(36) 



where 5* stands for the complex conjugate of 5. One of the main difficulties in 
writing down equations of motion for fields with a non-zero spin, and thus with 
more than one components, is to decouple the differential equation for a particular 
component from the remaining ones. This is not always possible but for the axially- 
symmetric Kerr-like black hole that we consider here it is. The situation actually 
is very similar to the pure four-dimensional case studied by Teukolsky 133 , where 
it was shown that at least the radiative components, that carry all the information 
for the propagating field, do indeed decouple 134 . Rearranging Eqs. (33)-(34) and 
using the identity 

(D — 2p — p*)(S - 2t) = 2t)(D - 2p) $x , (37) 

we obtain the following decoupled equation for <J> 

(D-2p- p*)(A + fj, - 2 7 ) $ - (6 ~ 2t)(5* + tt - 2a) $0 = . (38) 

Using the explicit forms of spin coefficients and operators, Eqs. (30) and (31), as 
well as the factorized ansatz (32), the above can be conveniently separated into an 
angular equation, 



f 



sin (9 d6 



d ( . dS?j 

— sm6> 



de 



2m cot e m 2 



sin e sin 2 9 



1 2 2 2 /i 

+ a uj cos e 



-2au> cos 9 + 1 - cot 2 e + A 



1.; 



S£(0)=O, (39) 



where Ay is a separation constant, and a radial equation, 



1 A 

A dr 



dR x 
dr 



K 2 - iKd r A 



+ Aiur + (A" - 2) - A 



i2i(r)=0, (40) 
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where we have defined: 



K = (r 2 + a 2 ) 10 — am , A s j = X s j + a 2 uu 2 — 2amu> . 



(41) 



Fermion Fields (s = 1/2). For a massless two-component spinor field, the Dirac 
equation can be written as: 

(6* - a + tt) xo = (D - p) xi , (42) 
(A + A1 -7) X o = (^ + /3-r)xi. (43) 
Performing a similar rearrangement as in the case of bosons and using the identity 
(5 - a* - t + tt*)(D - p) X i = (D - p - p*)(S + 0- T )xi, (44) 
we find that xi is decoupled leaving behind an equation for xo 

{D-p- p*)(A + ft - 7 ) xo - (6 - a* - r + n*)(S* - a + tt) X o = . (45) 

The latter can be explicitly written as a set of angular and radial equations, having 
the form 



1 d 



sin (9 d6 



dS 



smf 



1/2J 



dff 



+ 



m cot 9 m 2 
sin 9 sin 2 9 



2 2 in 

+ a co cos 9 



1 1 2 

— au>cos9+- — - cot 9 + Xi, 

2 4 2 " 



^y 2 j(<?)=0, (46) 



and 

1 d f A 3/2 dR i/2 



K 2 -iK(d r A)/2 n 1 . . »» > . 

A 2 2J 



y/A dr \ dr 
respectively, with the same definitions for K and A S j as before. 



R 1/2 (r) = 0,(47) 



Scalar Fields (s = 0). In the case of a scalar field, its equation of motion can be 
determined quite easily by evaluating the double covariant derivative g^ v D^D v act- 
ing on the field. As in the purely four-dimensional case 135 , the use of the factorized 
ansatz (32) leads, similarly to the previous cases, to the following pair of separated 
equations 



1 d 



sin 9 d9 



sin 9 ■ 



dS, 



0,1 



+ 



m 
sin 2 9 



+ a 2 io 2 cos^ 9 + \oe 



(48) 



d / dRo\ f K 2 \ , , 
where Y m e — e tmv S™ e (6) are now the spheroidal wave functions 131 . 



(49) 
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Combining all the above equations derived for bosons, fermions and scalar fields, 
we may now rewrite them in the form of a master equation valid for all types of 
fields. The radial equation then takes the form: 

A" s |- (A* +1 ^ ) + ^ 2 -^ A + Alsujr + s (A » _ 2) _ As ^ Rg (r) = o )(50) 



while the angular equation reads 



sin 9 — + 



sin0 d0 V d9 



2ms cot 9 m 2 2 2 

^ \- a u cos 1 

sin sin' 1 6> 



— 2aswcos# + s — s cot + A s 



S?j=0. (51) 



The latter equation is identical to the one derived by Teukolsky 133 in the case of a 
rotating, Kerr black hole. The radial equation differs by the extra factor s (A" — 2) 
due to the fact that for our metric tensor this combination is not zero, contrary to 
what happens in the case of the 4-dimensional Kerr, or Schwarzschild, metric. 

In order to solve the radial equation (50), we need to know the value of the 
constant term A s j defined in Eq. (41), or alternatively the value of the separation 
constant X s j that appears in the angular equation (51). Actually, the separability of 
radial and angular parts in the equation of motion of a field in a Kerr-like black hole 
background comes with a price: the separation constant is a complicated function 
of both the energy lj of the particle and the angular momentum parameter a of 
the black hole. An analytic form can be found, for any value of s, in the limit of 
auo — > 0, in which case we may write 136 > 137 > 138 

X sj = -a(a + 1) + E /r S KO fe = 3(3 + 1) - s{s + 1) - j^—j au + ... , (52) 

where terms higher than linear have been suppressed due to their complexity. It 
becomes therefore clear that any attempt to solve analytically the radial equation 
(50) must necessarily take place in the low-energy and low- momentum limit. 

Finally, we may simplify further the radial equation (50) by making the redefi- 
nition R s = A~ S P S . In that case, the A"-term disappears and we obtain: 

A ,rf / Al ,dP s \ {K 2 -isKd r A „. - \„,, „ 

A Tr [ A 17 ) + { A + AlSUT Asj ) Ps(r) = ° ' (53) 

where now A S j = A s j + 2s. 

Having derived the equation of motion for fields with arbitrary spin s propa- 
gating in the four-dimensional spacetime induced on the 3-brane, we can now move 
on to discuss the emission of Hawking radiation directly on the brane during the 
spin-down and Schwarzschild phase of a small, higher-dimensional black hole. 
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5.1. Emission of Hawking radiation during the spin-down phase 

The study of the emission of Hawking radiation on the brane during the spin-down 
phase has not been as thorough as one would hope. The only such attempt has been 
made by the authors of Rcf. 68, who independently performed the derivation of the 
master equation described earlier b . The master equation was solved analytically in 
the limit of small aui and for the particular case of only one extra compact dimension, 
i.e. n = 1. A well-known approximation method 136, 119 > 120 < 121 > 139-147 wag usec j 
according to which the asymptotic solutions at infinity and close to the horizon are 
found, and then matched in an intermediate zone in order to construct the complete 
analytic solution. The same method had been used earlier in Ref. 148 to compute 
the spectrum of Hawking radiation of scalar fields emitted both on the brane and in 
the bulk, and later, in Ref. 149, for the study of the emission of fermions and gauge 
bosons on the brane; both analyses were done for the Schwarzschild phase of the 
black hole, therefore, a description of the aforementioned approximation method 
will be given in the next subsection. 

In Ref. 68, the greybody factors, as well as the radiation spectra, for the emission 
of scalars, fermions and gauge bosons on the brane during the spin-down phase were 
derived for n = 1. For various values of the angular momentum parameter a, the 
radiation spectra were plotted as a function of the energy w of the emitted particle. 
The presented results reveal a suppression of the power spectrum for scalar fields, 
as the black hole angular momentum increases, both at low and high energies. On 
the other hand, in the radiation spectra for fermions and gauge bosons, we observe 
an enhancement at low energies that turns to a suppression at high energies. In all 
cases, the power emission curves, for all values of a and s, lie far below the ones 
produced by using the high-energy geometrical optics limit value of the greybody 
factors. 

However, these analytic results were derived in the limit of small auj, therefore, 
they are bound to break down at high energies. The high energy behaviour of the 
power spectra presented are thus not trustable and an exact (numerical) analysis 
is necessary to determine the complete spectrum. By comparing the results, for 
a = and n = 1, produced in Ref. 68 with the ones derived in the four-dimensional 
case 119 , we conclude that the power spectrum of all fields seems to be enhanced. 
However, for a > 0, the suppression caused by the black hole angular momentum 
complicates the picture and the comparison with the four-dimensional case is not 
clear. Similarly, no conclusions can be drawn on the behaviour of the radiation 
spectra as n increases. Apart from an exact analysis, valid at all energy regimes, a 
complete study of the dependence of the emission spectra on the number of extra di- 

b A typographical error appears in the published version of Ref. 68: the coefficient s that should 
have multiplied the A"-term [or n(n— f)/ir n_1 , in their notation] is constantly missing from their 
equations for scalars (A5), fermions (A23) and gauge bosons (A26), as well as from the master 
equation (A28). Nevertheless, the results presented in Ref. 68 are still robust since they were 
derived for the particular case of n = 1 for which this term vanishes trivially. 
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mensions, is clearly necessary. Both of these tasks have been successfully performed 
in the case of the Schwarzschild phase 150 , and the results are reviewed in the next 
subsection. 



5.2. Emission of Hawking radiation during the Schwarzschild phase 

In the four-dimensional case, the most optimistic estimates for the amount of the 
initial energy of the colliding particles lost in the balding (<~16%) and spin-down 
(~25%) phases leave behind approximately 60% of the total energy to be spent 
during emission of Hawking radiation in the Schwarzschild phase. Although we 
now know that these numbers change with the addition of extra dimensions, we 
believe that the Schwarzschild phase will still be the longer one and will account for 
the greatest proportion of the mass loss in the life of a small, higher-dimensional 
black hole. Contrary to the case of the spin-down phase, the emission of Hawking 
radiation on the brane during the Schwarzschild phase has been thoroughly studied, 
both analytically 148 - 149 and numerically 150 . In what follows we present a review 
of both types of results (see also 151 ). 

A higher-dimensional black hole with vanishing angular momentum induces on 
the 3-brane the following spherically-symmetric, Schwarzschild-likc line-clement 

ds 2 = -h(r) dt 2 + hir)- 1 dr 2 + r 2 (d9 2 + sin 2 dip 2 ) , (54) 

where still 

(\ n+l 
f ) ■ (55) 

The above line-element follows easily from the one in Eq. (22) after projecting out 
the extra angular coordinates, by fixing their values, and then setting a = 0. As 
in the case of the projection of a Kerr-like line-element on the brane, the induced 
metric still has an explicit dependence on the number of extra dimensions that exist 
in the fundamental, higher-dimensional theory 

The master equation (53) also takes now a much more simplified form, namely 



A s — (^A 1_s — ) + 
dr V dr 



r LU 2 r 2 „. isujr 2 h' 



+ 2iujsr A 



h h 



P s (r)=0 (56) 



where now A = A S j = j(j + 1) — s(s — 1). The above equation needs to be solved 
over the whole radial domain. We will first review the results that follow from the 
analytic approach 148 > 149 and then move on to discuss the ones following from the 
numerical analysis 150 . 

As mentioned in the previous subsection, the analytic treatment demands the 
use of an approximate method in which the above equation is solved at the near- 
horizon regime (r ~ th) and far-field regime (r 3> rjj), and the corresponding 
solutions are matched at an intermediate zone. In solving analytically Eq. (56), the 
spin s will remain an arbitrary parameter of the differential equation; this will allow 
us to solve the radial equation for particles with different spin s in a unified way. 
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Starting from the near-horizon (NH) regime, a change of variable r — > h(r) 
brings Eq. (56) to the form 149 



h(l-h) 



d 2 P s 
dh 2 



dP^ 
dh 



+ 



{ur H y 



+ 



2isojrn — A 



ISOJTh 



(n + l) 2 /i(l - A) (n + l) 2 (l-/i) (n+1) A 
In the above, we have made use of the relation 

dh _(n + 1) 



(l-h) 



P s = 0. (57) 



(58) 



that follows from the definition (55). Making a further redefinition of the radial 
function, P s (h) — h a (l — h)^F s {h), the above equation takes the form of a hyper- 
geometric equation for the radial function F Sl i.e. 



ab F, = . 



The hypergeometric indices (a, b, c) can be identified as 
s + n (1 — s) 



a = a + /3 + 



6 = a + /3, 



c = 1 - s + 2a . 



(59) 



(60) 



(n+1) ' 

The power coefficients a and /?, in turn, are found by solving second-order algebraic 
equations leading to the results 



a + — s + 



n+1 



n+1 



and 



Hi 



1 



1 - 2s ± ^{1 + 2j) 2 - 4Lu 2 r 2 H - 8isu>r H 



(61) 



(62) 



2(n+l) 

respectively. Then, by using the general solution 152 of the hypergeometric equation 
(59), together with the aforementioned relation between P s (h) and F s (h), we ob- 
tain the following general solution for the radial function P s (h) in the near-horizon 
regime: 

P NH (h) - A_h a (1 - hf F(a, b, c; h) 

+ A+ h~ a (1 - hf F(a - c + 1, b - c + 1, 2 - c; h) , (63) 

where A± are arbitrary constants. For simplicity, in the above, we have dropped 
the spin index s from the radial function P. 

An important boundary condition must be imposed on the general solution (63) 
at the horizon of the black hole: since nothing can escape from the black hole, in 
the limit r — > r H the above general solution must contain only incoming modes. 
In order to impose this condition, we expand Eq. (63) in the limit r — > ru, or 
equivalently h — ► 0, and we obtain 

P NH ~A^h a + A+ hr a = A- h s cxp(^r" +2 y) + A+ h- s cxp(-^r£+ 2 y) (64) 
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for a = a+, and 

P NH ^ A_ cxp(-iwr^ +2 y) + A+ cxp(wjr^ +2 y) (65) 

for a = a- . In the above, we have used the 'tortoise' coordinate y defined in Eq. (16). 
The aforementioned boundary condition at the horizon demands that A_ = 0, in 
the former case, and A + = 0, in the latter case. However, the remaining term in the 
asymptotic solution for a = a + then describes an incoming wave with a diverging 
amplitude at the horizon. This is clearly an irregular solution, therefore the choice 
a = a + must be discarded altogether. On the other hand, the choice a = a_ leads 
to a regular incoming wave with amplitude unity at the horizon. Finally, we need 
to make a choice for the sign appearing in the expression of the (3 coefficient. One 
can easily see that the criterion for the convergence of the hypergeometric function 
F(a, b, c; h), i.e. Re (c — a — b) > 0, clearly demands that we choose /? = /?_. 

Having found the solution in the near-horizon regime, we now turn our attention 
to the far-field asymptotic regime. In the limit r>rg,or/i^l, Eq. (56) assumes 
a simplified form, which reads 149 

d 2 P 2(1 - s) dP ( 2 2isuj A\ () 
dr 2 r dr \ r r 2 J 

Another redefinition of the radial function, P = e~ tur r J+;s P(r) 1 and a change of 
variable, z — 2iujr, puts the above equation in the form of a confluent hypergeo- 
metric equation 

d 2 P dP 

with a = j — s + 1 and b = 2j + 2. The general solution of the above differential 
equation is given by 152 

P{z) = B+M(a,b,z) + B-U(a,b,z), (68) 

where M and U are the Kummer functions, and B± are arbitrary coefficients. Then, 
the complete solution for the radial function P(r) at r > r/j can be written as 

P FF {r) = e-' lujr r j+s \b+ M(j -s + 1, 2j + 2, 2iujr) 

+B_ U(j -s + 1, 2j + 2, 2iLur)] . (69) 

In order to construct a complete solution for the radial function P(r), we need 
to match the two asymptotic solutions at an intermediate zone. Although they look 
different, we will see that they actually take a very similar form in the intermediate 
regime, which facilitates considerably their matching. In order to see that, we first 
shift the hypergeometric function towards large values of r. This can be done by 
using a standard linear transformation formula 149 > 152 ; that changes its argument 
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from h to 1 — h. Then, expanding in the limit r — ► oo, or h — > 1, we obtain 

ij+s r(l- S + 2a)r(-2/3+i n ^) 



iVJf 



(/') 



r « y r(a-/? + l-*)r(a-/3 + ±=^) 



+ ~^ r ' T(a + /3)r(a + l3+ s+n n %- a) 



(70) 



We then expand the far-field solution (69) in the limit urC 1, and we take 

Pff {?) = B + r*+> + -± T{2j + 1} (71) 

FF{ ' + + rJ- s +! r(j - s + 1) (2iw) 2 J+ 1 ' 1 J 

The matching of the two solutions (70) and (71) clearly provides relations between 
the integration constants appearing in the expressions of the asymptotic solutions. 
These read 



B, = 



A _ r(i-* + 2a)r(-2/? + ±=^) 



1 r(l - s + 2a) r(20 - i^) T(j -a + 1) 



B_ = 



A_ r° H s+1 (2ilu) 



T(a + (3) r(a + /? + ^+^1) T(2j + 1) 



(72) 



(73) 



The above relations complete the determination of the solution for the radial func- 
tion P s (r) that describes the propagation of scalars, fermions and gauge-bosons in 
the background of a projected (4 + n)-dimensional Schwarzschild-like black hole on 
a 4-dimensional brane. 

We may now proceed to the calculation of the greybody factors for each type 
of field. For this, we need the amplitudes of the incoming and outgoing modes at 
infinity. This can be found by taking the limit r — > oo of the far-field solution (69). 
Then, we obtain 



P^\r)=A 



(oo) 



(2ur) 



l-2s 



A 



(oo) 



+ 



with 



,(oo) 



-iw(j-s+l)/2 

{2uj)i+ s 



B- + 



° ut (2wr) 
B+ e^W-s+i) r(2j 



2) 



r(j 



i) 



,(oo) _ B + e-W s + 1 )/ 2 r(2j + 2) 



r(i 



1) (2c;) 



(74) 

(75) 
(76) 



c In order to obtain the simple powers of r shown in the expression below, the low-energy limit, 
ijjrn <C 1, has been taken in the expression of the (i coefficient in the first term of Eq. (63). No 
expansion has been made in the arguments of the Gamma functions. 
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A few comments are here in order: for s = 0, the above expression describes 
the incoming and outgoing mode of a scalar field at infinity, both scaling as 1/r; 
this expression, therefore, contains the complete information for the propagating 
scalar field and the two amplitudes, A^f 1 and A 1 ^ , can be used to determine the 
reflection coefficient 1Z, and from that the absorption coefficient A, as follows 

2 

(77) 

On the other hand, for s =/= and s = +\s\, Eq. (74) describes a dominant incoming 
mode and a suppressed outgoing one; vice versa, for s — —\s\, the outgoing mode 
dominates while the incoming one is greatly suppressed. This is merely a reflection 
of the fact that, for s ^ 0, the propagating field has more than one components, 
each one of which mainly describes cither the incoming mode or the outgoing mode. 
What is fortunate is that we do not need to know the complete solution at infinity 
as long as we can compute the total incoming flux Ti n at the horizon and at infinity. 
Then, the absorption coefficient may be directly determined through the following 
expression 

l-4 S) | 2 = ^7- (78) 

•'in 

The calculation of the incoming flux is indeed possible for the spherically-symmetric 
Schwarzschild-like background that we consider here d , so the only component whose 
expression we need to calculate, for fields with a non-zero spin, is the upper one, 
with s = +\s\. 

The incoming flux of a fermionic field can be computed from the radial com- 
ponent of the conserved current, = \[2 <7 AB ip A ip B , integrated over a two- 
dimensional sphere, first at infinity and then at the horizon 146 > 149 . For gauge bosons, 
the (ir)-componcnt of the energy momentum tensor T^ v = 2a AA ,a BB ,ip AE ''{p A B 
can be used instead. The same method for the calculation of the absorption coef- 
ficient can be equally well applied also in the case of a scalar field with the use of 
the radial component of the conserved current = hr 2 (tpd^tp* — ^*d^tp) - this 
method is completely equivalent to the one described above, that makes use of the 
expression of the reflection coefficient, while it allows us to write down a unified 
expression for the absorption probability for fields with arbitrary spin. Substituting 
the results for the incoming fluxes in Eq. (78), we may finally write 

|_4 S) | 2 = {2ur H f^ 



\A 



(0),2 



(0)|2 



A 



(oo) 



,(oo) 



IV"- 1 



.(oo) 



(79) 



d Note that for more complicated gravitational backgrounds, the calculation of the flux can 
be highly non-trivial, especially close to the horizon. In that case, the complete solution (both 
incoming and outgoing modes) at infinity must be determined. 
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where = A_, and A^f 1 is given in Eq. (75). In the above, we have taken into 
account the fact that, for fermions and gauge fields, the lower component contributes 
little to the in-falling flux both at infinity and at the horizon. We should note here 
that the coefficient in front of the ratio of the two amplitudes strongly depends on 
the field normalization (74) at infinity - although this is merely a convention, care 
should be taken so that this coefficient correctly reflects the chosen normalization e . 

Substituting the expression of A { ™ ] in Eq. (79) and using the relations (72) and 
(73), we obtain (for more information on the mathematical details omitted below, 
see Rcf. 149) 

UW I* = (2ur H )^ 

1 J 1 \T{l-s + 2a)\ 2 \C(ur H )^ +1 +D\ 2 ' 1 ' 

where the coefficients C and D stand for 

2 2j+i ^ (s-i/2) r f 2/3 - 1=£>) r(j - s + 1) 
C= t ^ ^ 

r(2. ? + 2)r(-2/3 + i^) 

D = ^ \l . (82) 

T{a-p + \-8)T\a-0+ ) T(j + s + l) 

Since we study the emission of brane-localized modes, this is clearly a four- 
dimensional process, therefore the relation (17) [or (18)] between the greybody 
factor and the absorption probability takes the simplified form 

-Sm - £ w + 1) i4 s) i 2 = (2i + 1) i4 s) i 2 . (83) 

where Ah = A-ktJj, or finally 

is) , , = (2ur H f^ (2j + l)A H 

j ' nK ' \T{1- s + 2a)\ 2 \C{ujr H f3+ 1 +D\ 2 ' 1 ' 

The coefficients C and D in the above expression depend, through Gamma func- 
tions, on the parameters of the theory, therefore the dependence of the greybody 
factor on the spin s, the energy w, or the number of extra dimensions n is not clear. 
Nevertheless, the above expression can be easily plotted and thus reveal the desired 
dependence. Before however discussing those results, we might attempt to derive 
some simple, elegant expressions for the greybody factor by taking the low-energy 
limit ojth -C 1. In this limit, we may express each one of the coefficients C and D 
as a power series in {loth) and keep only the leading term. Then, we may see 149 
that the term proportional to C, in the denominator of Eq. (84), is sub-dominant 
compared to the one proportional to D - that comes out to be independent of (corn) 

e For this reason, different coefficients appear in the expressions of the absorption probabilities 
1-4^ | 2 in Rcfs. 149 and 150, where different normalizations for the field at infinity were used. 
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at leading order - and thus the former can be ignored. Depending on the spin of the 
particle, additional powers of {uorn) may come from the measure of the remaining 
Gamma function in the denominator of Eq. (84). While for scalars and fermions, 
the zeroth-order term is (ur#)-independent, i.e. 

1 2ir 2 (iur H ) 2 



|r(l + 2a)| 2 3(n+l) 2 
1 1 2tt (luvh) 2 



+ OK) 4 , (85) 

lrd/2 + 2.)! 2 -. ' (n + 1) 2 + °^ 4 ' ™ 
for gauge bosons this is not the case and extra powers of {uorn) come up 

Putting everything together, we obtain the following simplified, low-energy expres- 
sions for the greybody factors 

(0) 7r(2j + i) r(^) 2 r(i + ^) 2 , (m 



for scalars, 



(89) 



(1/2) _ ,(2 J + l)2-(^ 2 )/(" +1 ) /WN*" 1 

,n [UJ) - 4r( J + i) 2 { 2 J Ah + -> 

fermions, 

(« + i) 2 ^ r ( 2 ^)r(2. ? + 2) y 

and gauge bosons. By choosing different values for j > s, the above approximate 
expressions reveal that, at the low-energy regime, the greybody factor for all types 
of fields gets suppressed as the angular momentum number increases; on the other 
hand, the same expressions predict, independently of the value of the spin s, an 
enhancement of the value of the greybody factor as the number of extra dimensions, 
that are projected on the brane, increases. 

Let us note at this point that the full analytic expression Eq. (84) has been 
derived in the low-energy approximation: the low-energy limit was taken in the 
expression of the f3 coefficient, during the matching of the two asymptotic solutions. 
The simplified analytic expressions, Eqs. (88)-(90), have been derived in the same 
low-energy limit by further expanding Eq. (84). However, it is clear that the more 
we expand an expression, the more we limit its validity, as more and more terms are 
left out. The simplified expressions (88)-(90), therefore, are valid in a more narrow 
low-energy regime than the full analytic one (84). Moreover, the information that 
we obtain from the simplified expressions is trustable only as long as it is deduced 
from the lowest partial wave, j = \s\, which always gives the main contribution to 
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Fig. 1. Analytical results for the greybody factors for emission of (a) scalars, (b) fermions, and 
(c) gauge bosons from a (4 + ra)-dimcnsional black hole on the brane. 



the greybody factor. The aforementioned dependence of cx^„(w) on j for all types 
of fields, as well as the dependence on n for fermions and gauge bosons were indeed 
derived by looking at the leading term of the lowest partial wave. However, in the 
case of scalar fields the dominant partial wave (j = 0) is found to be independent 
of n having the form a^(u)) ~ Attt^ — Ah, which is a well-known property of the 
low-energy behaviour of the greybody factor for scalar fields. One then is forced 
to look at higher-order terms for the desired dependence on n. Because the higher 
partial waves are enhanced with n, one may naively conclude that the greybody 
factor itself increases. However, it turns out that the next-to-leading order terms in 
the expansion of the lowest partial- wave, denoted by ellipses in Eq. (88), are of the 
same order as the leading terms of the higher partial waves, and thus should also 
be taken into account. When this is done, we are finally led to a decrease of the 
greybody factor for scalar fields with n. An extreme caution is therefore necessary 
when one deals with simplified expressions as the Eqs. (88)-(90). 

On the other hand, for more accurate, and thus more informative, evaluations of 
the greybody factors and associated emission rates, we clearly need to use the full 
analytic result (84). Although an approximate result itself, this expression leads to 
a behaviour for the greybody factors that is in excellent agreement with the exact 
behaviour in the low-energy regime, and in fairly good agreement with the one in the 
intermediate-energy regime. Plotting Eq. (84) as a function of the energy parameter 
ujrn, for different values of s, leads to the behaviour depicted in Figures l(a,b,c) 
149 , where a^ s is the greybody factor (84) summed over j up to the third partial 
wave - it is easy to see that any higher partial wave has a negligible contribution to 
the final result in this energy regime. As mentioned above, the greybody factor for 
scalar fields is found 149 to decrease with n, while the ones for fermions and gauge 
bosons are enhanced as n increases, at least up to intermediate energies. While 
scalars and fermions have a non- vanishing greybody factor as torn — * 0, the one for 
the gauge bosons vanishes. 

The full analytic result for the greybody factor (84) may then be substituted 
in the expression of the power flux for emission in the four-dimensional spacetime, 
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Fig. 2. Analytical results for the energy emission rates for (a) scalars, (b) fcrmions, and (c) 
gauge bosons from a (4 + n)-dimcnsional black hole on the brane. 



namely 



-^T~ ^ ^> (W) exp(u,/T ff )±l 2^ ' (91) 



Figures 2(a,b,c) depict the behaviour of the energy emission rates for particles with 
spin 0, \ and 1 in the low- and intermediate-energy regime. Despite the differ- 
ent low-energy behaviour of the greybody factors for scalars, fermions and gauge 
bosons, depicted in Figure 1, the corresponding power emission rates, as well as the 
flux emission rates, exhibit a universal behaviour according to which the energy, 
and the number of particles, emitted per unit time and energy interval is strongly 
enhanced, as n increases. This was to be anticipated: we remind the reader that 
the temperature of the black hole is given by the relation Tjj = (n + l)/47rr£r, 
therefore, for fixed r#, the temperature of the black hole increases as n increases. 
This simply means that the energy of the black hole available for the emission of 
particles also increases, and this is reflected in the enhancement of the power and 
flux rates. The effect of the greybody factor, on the other hand, is to suppress or 
enhance - depending on the spin of the particle, the dimensionality of spacetime 
and the energy regime that we are looking at - the emission rates compared to the 
ones derived by using its high-energy, geometrical optics value. This suppression or 
enhancement is then reflected at the location of the peak of the emission curve. For 
example, having seen that the greybody factor for gauge bosons is very much sup- 
pressed at the low-energy regime, compared to scalars and fcrmions, we expect that 
more spin-1 particles will be emitted in the intermediate and high-energy regime, 
thus shifting the peak of the gauge boson emission curve towards higher energies. 
This feature will indeed be obvious when the complete spectra are derived. 

Each one of the curves depicted in Figure 2 has a blackbody profile. Neverthe- 
less, the approximate analytic expression (84) for the greybody factor, and thus 
the corresponding emission rates, are accurate only in the low-energy regime while 
at intermediate energies it mainly provides qualitative agreement with the exact 
behaviour. Any attempt to derive emission rates in the high energy regime by using 
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this expression is bound to lead to wrong results. For this reason, the graphs in 
Figure 2 extend only up to intermediate energies for fermions and gauge bosons, 
while for scalars, our analytic formula breaks down much earlier [see Fig. 1(a)]. 

Since the analytic results can accurately describe only the part of the "greybody" 
curve that extends over the low-energy regime, a numerical analysis is clearly nec- 
essary in order to derive the complete radiation spectrum. In that case, the master 
radial equation (53) can be solved through numerical integration 150 under the same 
boundary conditions that were used in the analytic approach: a purely incoming 
mode at the horizon of the black hole is integrated outwards until the asymptotic 
solution (69) is reached. The determination of the asymptotic coefficients, A_ and 
, can then be used in conjunction with Eq. (79) to define the absorption co- 
efficient, and from that the greybody factors and emission rates 150 , as described 
above. The numerical manipulation of this problem is not without technical difficul- 
ties itself: as the spin of the particle increases, it becomes more and more difficult to 
accurately determine the two asymptotic coefficients; employing various transfor- 
mations of the radial equation can help remedy this problem (for more information 
on this point, the reader is referred to Ref. 150). 

The numerical integration was performed in Ref. 150 where both greybody fac- 
tors and emission rates for particles with spin 0, \ and 1 were computed. The exact 
results for the greybody factors are displayed in Figures 3(a,b,c). Their complex 
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Fig. 3. Numerical results for the greybody factors for emission of (a) scalars, (b) fermions, and 
(c) gauge bosons from a (4 + n)-dimcnsional black hole on the brane. 
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behaviour depending on the spin of the particle emitted, the number of the extra 
dimensions projected on the brane and the energy regime studied, is evident. At the 
low-energy regime, the exact results closely follow the ones derived by the analytic 
method. As loth — > 0, the greybody factor for scalar fields reduces to Airrfj, thus 
revealing the fact that the low-energy value of this quantity is indeed proportional 
to the area of the black hole horizon, even when a number of extra dimensions is 
projected onto our four-dimensional spacctime. The greybody factor for fermions 
also reduces to a constant value, different for different values of n, while the one for 
gauge bosons goes to zero as predicted by the analytic results. 

As the energy increases further, the exact results start deviating from the ana- 
lytic ones and they can be seen adopting an oscillatory behaviour - due to the late 
dominance of higher partial waves - around a high-energy asymptotic value. This 
asymptotic value is the same for all particle species, nevertheless it strongly depends 
on the number of extra dimensions. In this regime, the greybody factor assumes its 
geometrical optics limit value, already well-known from the four-dimensional case 
119, 121,123, 124_ -p Qr a mass i ess particle in a circular orbit around a black hole, de- 
scribed by the line-element (7), its equation of motion, p^p^ = 0, takes the form 

rdtp) & 2 r 2 ' [ ' 

where b is the ratio of the angular momentum of the particle over its linear momen- 
tum. The classically accessible regime is defined by the relation b < min {r/y/h). Al- 
though the number of the projected dimensions do not change the general structure 
of the above equation, valid for motion in a four-dimensional, spherically-symmetric 
background, the expression of the metric function h(r) does change as it carries an 
explicit dependence on n. Using the definition (55), we can easily find that the 
closest the particle can get to the black hole is at a distance 155 



l/n+l 



The above expression correctly reproduces, for n = 0, the four-dimensional value 
of r c = 3\/3?"£r/2 121,123 . The radius r c defines the absorptive area of the black 
hole at high energies and, thus, the corresponding value of the greybody factor, 
a g = 7rr 2 , which, being a constant, now describes a blackbody. Since r c decreases 
as n increases, the asymptotic greybody factor becomes more and more suppressed 
as the number of extra dimensions projected onto the brane gets larger. 

Turning to the emission rates, the exact results for the energy emission in the 
form of scalars, fermions and gauge bosons are depicted in Figures 4(a,b,c) 15 °. The 
complete emission curves have now been constructed and can be easily compared for 
different values of n and for different species. For all types of particles, an important 
enhancement of the energy emission rate takes place, as n increases. This was already 
predicted by the analytic results, however, the exact numerical results, that are now 
available, allow us to accurately calculate the enhancement of the emission rates. 
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Fig. 4. Numerical results for the energy emission rates for (a) scalars, (b) fcrmions, and (c) 
gauge bosons from a (4 + n)-dimensional black hole on the brane. 

Integrating over the whole energy regime, we obtain the total energy emissivities, 
which arc displayed in Table 4 150 normalized to the corresponding emissivities for 
n = 0. It is worth pointing out that the amount of energy emitted by the black hole 
in the case of n = 7 is 3 orders of magnitude larger than the one for n = 0, for 
all types of particles. A similar enhancement is observed in the number of particles 
emitted by the black hole. 



Table 4. Total energy emissivities for different values of n 



n 





1 


2 


3 


4 


5 


6 


7 


Scalars 


1.0 


8.94 


36.0 


99.8 


222 


429 


749 


1220 


Fcrmions 


1.0 


14.2 


59.5 


162 


352 


664 


1140 


1830 


Gauge Bosons 


1.0 


27.1 


144 


441 


1020 


2000 


3530 


5740 



Another important question, that one can pose, is what type of particles the 
black hole 'prefers' to emit, as the number of extra dimensions n changes. From the 
four-dimensional analyses 119 ' 124 we know that, for n = 0, most of the energy of the 
black hole is emitted in the form of scalar particles; this becomes clear by comparing 
the areas under the emission curves, for particles with spin 0, \ and 1, in Figure 
5(a). The emission curves in Figure 5(b) have been drawn for n = 2: the emission 
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Fig. 5. Energy emission rates for scalars, fermions and gauge bosons for (a) n = 0, (b) n = 2, 
and (c) n = 6. 

curves for scalars, fermions and gauge bosons now seem to define comparable areas. 
Increasing further the value of n, i.e. n = 6, we obtain the emission curves shown 
in Figure 5(c). The situation has radically changed: the gauge bosons are now the 
dominant 'channel' for the emission of the black hole energy, with the scalars and 
fermions coming second and third, respectively. Accurate estimates for the relative 
emissivities have been made 150 , and are displayed in Table 5, normalized to the 
emissivity of scalar fields. The 3-species ratios shown below can equally well be 
used as signatures for the dimensionality of spacetime if the emitted spectrum of 
Hawking radiation is successfully observed. 



Table 5. Relative energy emissivities for different values of n 



n 





1 


2 


3 


4 


5 


6 


7 


Scalars 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


Fermions 


0.55 


0.87 


0.91 


0.89 


0.87 


0.85 


0.84 


0.82 


Gauge Bosons 


0.23 


0.69 


0.91 


1.00 


1.04 


1.06 


1.06 


1.07 



6. Hawking Radiation in the Bulk 

We now turn to the study of the Hawking radiation emitted by a small, higher- 
dimensional black hole in the bulk. According to the assumptions of the model, 
only gravitons and scalar fields can propagate in the space transverse to the brane, 
therefore, these are the only modes that should be considered here. Although this 
type of emission is unlikely to be observed, the study of the bulk emission is ex- 
tremely important since it determines the amount of energy lost in the bulk, and 
subsequently the amount of energy left for emission in the form of brane-localizcd, 
observable modes. 

Due to the simplicity of the metric, the emission of modes from a spherically- 
symmetric, higher-dimensional black hole was studied first 148 . In addition, the 
absence of any technical difficulties in the derivation of the equation of motion 
for scalar fields propagating in a higher-dimensional spacetime naturally led to the 
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study of this particular type of modes first. The emission of bulk scalar modes 
by a Schwarzschild-like black hole has been studied both analytically 148 ' 153 and 
numerically 150 . We start by reviewing these results before presenting the definite 
answer to the question of bulk-to-brane relative emissivity raised above, and before 
commenting further on the emission of scalar and graviton modes in the bulk. 

The gravitational background around a (4 + n)-dimensional, spherically- 
symmetric black hole is described by the line-element (7) 113 , which we give again 
here for convenience 



ds 2 



-h(r) dt 2 + h(r) 



r 2 dn 2 +n . 



(94) 



The equation of motion of a scalar field propagating in the above background is 
simply given by DmD m cf> = 0, where Dm is the covariant derivative in the higher- 
dimensional spacetime. This equation takes a separable form if we use the ansatz 



<f>(t,r,0i,<p) = e- vat R ue (r)Y i (Sl), 



(95) 



where Yi{Q) is the (3 + n)-spatial-dimensional generalization of the usual spherical 
harmonic functions depending on the angular coordinates 116 . The radial equation 
then takes the form 148 



hjr) d 

r n+2 ^ r 



h{r) r 



n+2 



dR 

dr 



+ 



R = 0. 



(96) 



As in the case of emission on the brane, the above equation must be solved and 
the absorption coefficient must be determined before we are able to write down 
the expression for the greybody factor and the emission rates. The same analytic 
method was also used in this case 148 : Eq. (96) was solved in the near- horizon and 
far-field regime and the two solutions were matched in an intermediate zone. In 
the near- horizon regime, the same change of variable, r —* h(r), brings the radial 
equation to the form 

{ur H ) 2 e(£+l + n) 



i j n d 2 R dR 



R = 0, (97) 



(n + l) 2 h(l - h) (n + 1) 2 (1 -h) 

while the redefinition R{h) = h a (l — hYF{h) reduces it to a hypergeometric equa- 
tion with a = b = a + [3 and c = 1 + 2a, where 



a± = ± 



(98) 



n + 1 ' 

Then, the general solution of Eq. (97) is written as 152 

R NH (h) = A_h a± (1 - hf F(a, b, c; h) 

+ A+ h- a± (1 - hf F(a - c + 1, b - c + 1, 2 - c; h) . (99) 

If we expand the above solution in the limit r — > r#, or h — > 0, and choose a = a_, 
we obtain the approximate form 



R NH ~ A_ exp(-^r" +2 y) + A+ cxp(iujr^ +2 y) , 



(100) 
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in terms of the tortoise coordinate (16). Had we chosen the alternative option 
a = a + , the same form would have been derived with the coefficients A_ and 
A + interchanged. Since the two options are completely equivalent, we choose the 
former one, and we demand A + = in order to have a purely incoming mode near 
the horizon of the black hole. The criterion of the convergence of the hypergeometric 
function Re (c — a — b) > demands once again that /? = /?_. 

In the far-field zone, the limit r^>rn and the redefinition R(r) — f{r)/A n+1 ^/ 2 
reduce Eq. (96) to a Bcssel differential equation of the form 148 

d 2 f 1 df 



dr 2 r dr~^ 



u 2 



-L(, + H±i) 2 ] / = . (ioi) 



152 



The general solution then for the radial function R(r) is given by 

B B 
R Mr) = r(n+ t)/2 J £+(n+i)/2(ur) + r(n+1)/2 Y e+(n+1)/2 (ujr) , (102) 

where J and Y are the Bessel functions of the first and second kind, respectively. 

In order to match the two asymptotic solutions, the far- field one is expanded in 
the limit ur <C 1, giving 

n B+r e S_ / 2 \ e+{n+1)/2 r(^+2±l) 

while the near-horizon solution, after being "shifted" and its argument changed to 
1 — h 152 , is expanded in the limit r » r# thus assuming the form 

'fry r(l-2/3) /r H y+«+i T(2/3- 1) 
_feJ r(l + a-/3) 2 + \~) T(a + f3) 2 

Matching the two solutions (103) and (104), we obtain a relation between the two 
integration constants at infinity 

o / 2 n 2£+n+i r(f + 2±i)2 + «±i) r(i - 2/?) r(« + /3) 2 

' > ^ Z . (105) 



R NH {h)~A_T{\ + 2a) 



(104) 



B_ V^r^/ 7rr(l + a-/3) 2 r(2/3- 1) 

After having completed the determination of the solution for the radial function 
R(r), we turn our attention to the form of the scalar field at infinity. We need to 
determine the amplitudes of the incoming and outgoing modes, thus, we expand 
Eq. (102) in the limit r — > oo, and we find: 



with 



^(oo) = A oo) + A (oo) _J_ (106) 



= B ++^- e -(*+i +D/2 ; (10?) 

V2W7T 

4^ = e -^+t + D/ 2 . (108) 
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The reflection coefficient TZe is defined as the ratio of the outgoing amplitude over 
the incoming one at infinity. Then, the absorption coefficient At is given by 



|^| 2 = 1-|^| 2 = 1- 



B+ - iB. 



2i ^-^ (109) 



BB* + i(B* - B) + 1 



B+ + iB- 

whcrc B = B+/B- is defined in Eq. (105). 

The above analytic result can take a simplified form in the low-energy limit 
ujr H <C 1, in which case BB* ^> i(B* — B) » 1, and we may write 

, 2 



/ x2t+n+2 rfl + — — Y 



^ ,110 ' 



Substituting the above result into Eq. (18), we arrive at 

2 



148 



/ \2i r(i + — —\ r(^-\ 

( I, 1 + n+1 J { 2 J 



»'.»(") = SZTKIiT I — ) -T 1 ^ ««/!„ + ..., (Ill) 

where is the multiplicity of states with the same angular momentum number t 
defined in Eq. (19). Since loth <C 1, the greybody factor decreases as t increases, 
therefore, the main contribution to ai >n comes from the lowest partial wave 
with t = 0. It is easy to see that the above expression evaluated for I = simply 
reduces to Ah, thus, revealing the fact that even in the higher-dimensional case, 
the greybody factor for scalar fields at the low-energy regime is given by the area of 
the horizon f . This behaviour is similar to the one obtained in the four-dimensional 
case; here, however, the area of the horizon changes as n varies. 

As in the four-dimensional case, the contribution to the greybody factor from 
the dominant partial wave comes out to be independent of the number of extra 
dimensions. Looking at the dependence of the higher partial waves on n, we obtain 
a suppression of the greybody factor as the dimensionality of the bulk increases. 
However, in order to be absolutely certain about this behaviour we would have 
to include next-to-leading-order corrections in the simplified expression of ag tn (uj) 
(111), or simply deal with the full analytic result derived from Eq. (109). In either 
case, however, the derived dependence would only hold in the low-energy regime 
and no information could be derived from these expressions for the dependence of 
the greybody factor, and thus of the emission rates, in the high-energy regime. 

For this reason, we turn to the numerical integration of the radial equation (96) 
that can provide exact results for the greybody factor and energy emission rate in 
the bulk valid at all energy regimes. The behaviour of both these quantities are 
shown in Figures 6(a,b) 150 . As it was predicted by the analytic expression, the 
greybody factors, normalized to the area of the horizon, tend to unity for all values 

f This has been recently shown 154 to hold also for massive scalar particles propagating in a 
higher-dimensional spacctimc. 
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Fig. 6. (a) Greybody factors and (b) energy emission rates for scalar emission in the bulk from 
a (4 + n)-dimcnsional black hole. 

of n, at the low-energy regime. The naive prediction, also coming from the analytic 
expression, for the suppression of the greybody factor with n, at the same energy 
regime, is verified. Similarly to the emission on the brane, as the energy parameter 
ujrff increases, the greybody factor oscillates around an asymptotic high-energy 
limiting value. In analogy with the four-dimensional case, we expect the greybody 
factor to adopt the (4 + n)-dimensional geometrical optics limit at the high-energy 
regime. For large values of the energy of the scattered particle, the greybody factor 
becomes equal to the area of an absorptive body of radius r c which is projected on 
a plane parallel to the one of the orbit of the moving particle 123 . We may compute 
this area by setting one of the azimuthal angles 8i equal to 7r/2 and integrating over 
the remaining angular coordinates. Then, we find 150 

Ap ~ (^+2) T[(n + 2)/2] Tc ■ (112) 

The area of the absorptive body is clearly n-dependent, and reduces, as expected, 
to the usual four-dimensional result A p = 7rr 2 , for n = 0. In the above, r c is the 
value of the effective radius defined in Eq. (93) 155 . The greybody factor then is 
given by the above projected area, or more explicitly, 

1 r[(n + 3)/2] ^n + 3V" +2)/( " +1) /n + 3\ ( " +2)/2 
CTfl -^(« + 2) r[(n + 2)/2] {—) AH - 

In Ref. 150, it was pointed out that the high-energy asymptotic value of the grey- 
body factor, determined through numerical integration, although very close to the 
ones following from the above analytic expression, did not exactly coincide. How- 
ever, further study of the asymptotic behaviour revealed that the numerical results 
do indeed agree with the analytic prediction. The apparent disagreement was due to 
the fact that the bulk greybody factor adopts its asymptotic value at much higher 
energies than the brane greybody factor - extending our numerical calculation to 
higher energies led to the complete agreement with Eq. (113). 
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According to Figure 6(b), the emission rate of scalar fields in the bulk is en- 
hanced as the number of extra dimensions increases. As with the emission on the 
brane, this is caused by the increase in the temperature of the black hole, which 
eventually overcomes the decrease in the value of the greybody factor and causes 
the enhancement of the emission rate with n at high energies. We can now provide 
an accurate answer to the question of how much energy is emitted into the bulk 
compared to the one on the brane. By using Eq. (21), the relative bulk-to-brane 
energy emission rate may be written as 

dE^/dt TM B) \A B) ? nu) 

dsw/dt Y.M b) \^\ 2 ' 

where the superscripts B and b denote emission in the bulk and brane, respectively. 
In the above, we have taken into account that bulk and brane modes 'feci' the 
same temperature. The multiplicity of states in the bulk rapidly increases with n, 
compared to the one on the brane that remains constant. However, it turns out that 
the enhancement of the absorption probability for brane emission, as n increases, 
is considerably greater than the one for bulk emission. The exact numerical results 
derived in Ref. 150 allow us to compute, once again, the total emissivity of the black 
hole in the bulk and compare it with the one on the brane. The total bulk-to-brane 
relative emissivities, as a function of n, are given in Table 6 150 . 

Table 6. Bulk-to-brane energy emissivities for different values of n 



n 





1 2 


3 


4 


5 


6 


7 


Bulk/Branc 


1.0 


0.40 0.24 


0.22 


0.24 


0.33 


0.52 


0.93 



In Ref. 155, it was shown that the whole tower of KK excitations of a given 
particle carries approximately the same amount of energy as the massless zero- 
mode particle emitted on the brane. Combining this result with the fact that many 
more types of particles live on the brane than in the bulk, it was concluded that 
most of the energy of the black hole goes into brane modes. The results obtained 
in Ref. 155 were only approximate since the dependence of the greybody factor 
on the energy of the emitted particle was ignored and the (low-energy only valid) 
geometric expression for the area of the horizon was used instead s . The numerical 
results reviewed above have used the exact dependence of the greybody factor on 
both the energy and number of extra dimensions in order to provide the most 
accurate estimate for the bulk-to-brane energy emissivity. From the entries of the 
above Table, it becomes clear that the emission of bulk modes for intermediate 
values of the number of extra dimensions, i.e. n = 2,3,4 and 5, is particularly 
suppressed compared to the one of brane modes. For lower or higher values of n, 

g In Ref. 156, it was correctly argued that the dependence of the greybody factor on n must be 
taken into account, however, the dependence on energy u> was again ignored. 
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the emission of bulk modes becomes important but the bulk-to-brane ratio never 
exceeds unity. In conclusion, most of the energy of the higher-dimensional black 
hole, in the 'scalar' channel, is indeed emitted directly on the brane, nevertheless, 
depending on the number of extra dimensions, a substantial amount of the total 
energy may be lost in the bulk. 

However, a definite conclusion regarding the total amount of energy which is lost 
in the bulk cannot be drawn before the emission of gravitational radiation is also 
studied. A master equation describing scalar, vector and tensor gravitational per- 
turbations, in a higher-dimensional spherically-symmetric spacetime, was recently 
derived 157 and was shown to take the form of a 2nd order differential wave equa- 
tion. At the moment, the interest has been focused on the calculation 158 459,i60,66 
of the corresponding quasinormal modes 161 , resonances at complex frequencies 
that dominate at late times after the perturbation of the black hole background 
(see also 162 > 16 3a64,i65,i66 f or study of quasinormal modes of scalar fields in 
the higher-dimensional spacetime) . No study addressing the question of the relative 
bulk-to-brane emission rate for gravitons has yet been done. 

The emission of Hawking radiation in the bulk during the spin-down phase has 
again not been studied as much as the one during the Schwarzschild phase. The 
equation of motion of a scalar field in a five-dimensional Kerr-like background was 
shown 167 to take a separable form (just like the Hamilton-Jacobi equations for 
particles and light propagating in the same background 168 ) and analytic formu- 
lae for the energy and angular momentum fluxes, valid at the low-energy regime, 
were written down; no quantitative results, however, for the radiation spectra were 
produced and no generalization of this analysis for arbitrary number of extra di- 
mensions has been carried out. The study of the stability of higher-dimensional, 
rotating black holes has also attracted some attention: in Ref. 169 it was argued 
that 'ultra-spinning' (4 + n)-dimensional black holes, i.e. with an arbitrarily large 
angular momentum, are unstable - an earlier work had suggested that black holes 
with an angular momentum beyond a critical value will 'decay' to a rotating black 
ring 170 . On the other hand, perturbations on a massless scalar field propagating 
in a Kerr-like, higher-dimensional black hole was shown 171 to be free of unstable 
modes, no matter how large the angular parameter a is. Finally, the stability of a 
five-dimensional rotating black hole projected onto the brane was checked 172 , by 
studying scalar, electromagnetic and gravitational perturbations, and no unstable 
modes were found in the spectra. It is worth mentioning at this point that the 
Schwarzschild-likc higher-dimensional black holes were shown 173 to be stable un- 
der vector and tensor gravitational perturbations, with a potential instability arising 
only from the scalar gravitational sector, for n > 2. 

7. Conclusions 

The proposal of the existence of extra dimensions in the Universe has opened many 
different pathways that, in principle, lead to important modifications in the four- 
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dimensional cosmology, particle physics phenomenology and black hole physics. The 
recent revival of the idea of the existence of extra space-like dimensions, that can 
have an almost macroscopic size or be even non-compact, introduced a novel feature 
in the theory: the scale at which the gravity becomes strong may be much lower 
than the traditional four-dimensional Planck mass Mp. This has led to the exciting 
prospect that high-energy collisions between elementary particles, that will take 
place at next-generation ground-based accelerators or have been already taken place 
at the atmosphere of the Earth, can probe the energy regime of quantum gravity. 
Some of the most striking consequences would be the possible creation of mini black 
holes, or even strings and D-branes, as the products of a high-energy collision of 
particles with center-of-mass energies at just a few times the new scale of gravity 

In Section 2, we have reviewed the existing results in the literature regarding the 
possibility of the creation of mini black holes during high-energy collisions. Well- 
known, four-dimensional analyses have been generalized to cover the case where 
the colliding particles propagate in a higher-dimensional spacetime. Some of the 
new studies have shown that, in the case of head-on collisions, the mass of the 
produced black hole gets suppressed, as the number of extra dimensions increases, 
while some recent complimentary results indicate that the emission of gravitational 
waves during the collision is actually suppressed when n increases: this obviously 
brings the two types of results in an apparent disagreement unless we accept the 
possibility that a significant amount of the energy, lost during the head-on collision, 
is emitted in a form different from that of gravitational radiation. On the other 
hand, for collisions with a non-vanishing impact parameter, which are the most 
likely to take place, the black hole production cross-section is in fact enhanced with 
the number of extra dimensions. Putting the above into the framework of a realistic 
collision between composite particles have led to large estimates, by new physics 
standards, for the corresponding black hole production cross section, either at the 
LHC or at the atmosphere of the Earth. 

The effect of extra dimensions is not restricted to the production of mini black 
holes; the properties of the produced, higher-dimensional black holes are also mod- 
ified. In Section 3, we have discussed some of those properties, namely, the horizon 
radius, temperature and lifetime, all of which arc crucial parameters for the suc- 
cessful production and detection of these elusive, up to now, objects. As we saw, 
the horizon radius of a (4 + n)-dimensional black hole is many orders of magnitude 
larger than the one of a four-dimensional black hole with the same mass, which 
simply means that a mass M needs to be compacted less in a higher-dimensional 
spacetime to create a black hole. The temperature of these small black holes, in 
turn, comes out to be lower than in four-dimensions, which means that the emis- 
sion rate of Hawking radiation is smaller and their lifetime longer. What is most 
favourable for the possibility of detecting these objects is the fact that, for a black 
hole with a mass Mbh = 5 TeV or so, the Hawking radiation spectrum reaches its 
peak at energies close to the temperature of the black hole which is of the order 
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of 100-600 GeV (for n = 1, ...,7); this is exactly the energy regime that present 
and next-generation collider experiments can probe. Such a small black hole will be 
extremely short-lived, i.e. r ~ 10 -26 sec, nevertheless, the proximity of the evap- 
orating black holes to our detectors increases significantly the possibility of their 
detection. 

These mini black holes, upon implementation of quantum effects, emit Hawking 
radiation into the higher-dimensional spacetime in the form of both bulk and brane 
modes. The generalization of the four-dimensional expressions for the emission rates 
in the case of a higher-dimensional spacetime is straightforward, nevertheless, the 
exact expression of the greybody factors for different types of fields propagating 
in a higher-dimensional background was, until recently, unknown. As we explained 
in Section 4, the greybody factors encode valuable information for the background 
around the emitting black hole and depend on the energy of the emitted particle, 
its spin and the dimensionality of spacetime. This means that the presence of the 
greybody factor in the radiation spectrum will cause the modification of the low- 
energy emission rate from the high-energy one, and will lead to different emissivities 
for particles with different spin. Moreover, both the number and the type of particles 
emitted will depend on the number of extra dimensions that exist in nature, a feature 
that may possibly lead to the determination of the dimensionality of spacetime upon 
detection of Hawking radiation. 

Having concluded that the implementation of the greybody factor in the radia- 
tion spectrum is imperative in order to derive accurate estimates for the emission 
spectrum of the black hole, we moved on, in Section 5, to derive a master equa- 
tion for the propagation of fields with different spin on the black-hole background 
induced on the four-dimensional brane. By solving this master equation, one can 
compute the transmission cross-section, in other words the greybody factor, for 
brane-localizcd modes emitted by the black hole. This type of emission during the 
spin-down phase of a black hole has been only partially studied: radiation spectra 
for fields with different spin s have been computed only for the case of n = 1 in the 
limit of low energy and low angular momentum. On the contrary, the spherically- 
symmetric Schwarzschild phase has been thoroughly investigated. Both analytical 
and numerical methods were used, with the former one leading to analytical, but 
low-energy-only-valid, expressions, and the latter one yielding exact numerical re- 
sults valid at all energy regimes. We were thus able to compare the radiation spectra 
for different types of fields and different number of extra dimensions. For all types 
of particles, the total emissivities are greatly enhanced with the number of extra 
dimensions, with the enhancement reaching orders of magnitude for large values of 
n. As the increase in the emission rate depends strongly on the spin, the relative 
emissivities for particles with different spin are also strongly n-dependent: while 
scalar fields remain the preferred type of particle emitted by the black hole for low 
and intermediate values of n, they arc outnumbered by the gauge bosons for large 
values of n, with the fermions being the least effective channel during the emission. 
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The emission of brane-localized modes is undoubtly the most phenomenologi- 
cally interesting effect since it involves Standard Model particles that can be easily 
detected during experiments. On the other hand, the emission of bulk modes will be 
only perceived as a missing energy signal by the observer on the brane. Nevertheless, 
the amount of energy lost in the bulk is crucially important as it determines the 
remaining available energy for emission on the brane. The details of the emission 
of bulk scalar modes during the spin-down phase have been studied in an analytic 
but qualitative way, and no quantitative results are available. For the Schwarzschild 
phase, as we saw in section 6, the study has been completed. The greybody fac- 
tor and emission rates have been calculated and the latter were shown to increase 
again with the dimensionality of spacetime. The extremely important question of 
the bulk-to-brane energy emissivity has been answered only for the scalar channel 
and for the Schwarzschild phase: the amount of energy spent by the black hole for 
the emission of bulk scalar modes is always smaller than the one for brane modes; 
nevertheless, the amount of energy lost in the bulk must always be taken into ac- 
count especially for large values of n when the energies spent in the brane and 
bulk channel become comparable. No results for the bulk-to-brane emissivity for 
gravitons have been yet derived. 

Although the possibility of the production and evaporation of mini black holes 
at the LHC is an exciting prospect, this will be possible only in the case where the 
fundamental scale of gravity M* is indeed very close to 1 TeV. Nevertheless, there is 
absolutely no guarantee for that, and the only argument in favour of this particular 
value is the possible resolution of the hierarchy problem. If M* is larger than 1 TeV, 
even by one order of magnitude, the probability of the production of black holes at 
the LHC vanishes (although we might still witness these type of effects in cosmic 
ray particles). Nevertheless, all the analytical and numerical results presented in 
this review have M* as an independent parameter, and are therefore valid for all 
values of M*. If this scale is pushed upwards by various constraints, the derived 
results will still be applicable for the production and evaporation of black holes at 
the new energy regime \fs > M* . 

We would also like to stress that the results for the radiation spectra reviewed 
here refer to individual degrees of freedom and not to elementary particles, like 
electrons or quarks, which contain more than one polarization. For the number of 
elementary particles produced, and the energy they carry, one has to use a Black 
Hole Event Generator 71 > 174 . This simulates both the production and decay of small 
black holes at hadronic colliders and provides estimates for the number and spectra 
of the different types of elementary particles produced. 

In this review, we concentrated on theories postulating the existence of Large 
Extra Dimensions, and we studied the properties of small black holes that live in a 
spacetime with D — 1 fiat spacelike dimensions. Our analysis is definitely not valid 
in highly curved spacetimes, like the ones in five-dimensional warped models 3 . 
In that case, attempts to construct a gravitational background that would reduce 
to a black-hole line-element on the brane while being well defined away from it 
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have failed; numerical methods have instead found five-dimensional localized black 
holes that do not necessarily have a black hole line-element projection onto the 
brane (see Refs. 175-208 for some relevant works). Nevertheless, in the case of very 
small bulk cosmological constant, the 'warping' of the five-dimensional spacetime, 
parametrised by the inverse AdS radius, A -1 , is small and the extra spacetime can 
be considered as an almost flat one. Alternatively, if the horizon radius r# is much 
smaller than the AdS radius, no matter what the value of A is, then the black hole 
cannot perceive the warping of the fifth dimension. In either case, we may model 
black holes with m <C A arising in a warped spacetime as black holes living in a 
five-dimensional, flat spacetime. Under this assumption, all the results presented in 
this review hold equally well also for the 'warped' black holes. 

As a concluding remark, we would like to stress once again that the detection 
of signatures of possible black hole production events during high-energy collisions 
would be a revolutionary development both for particle physics phenomenology 
and gravitational physics. Any observational signal of this type would automati- 
cally prove the existence of extra dimensions and of a fundamental theory of all 
forces, with a low energy scale, valid in a higher-dimensional spacetime. The detec- 
tion of Hawking radiation emitted by these small, higher-dimensional black holes is 
the most direct evidence for the production and evaporation of those objects. The 
radiation spectrum of a decaying black hole can also reveal the exact dimensionality 
of spacetime as both the amount and type of the emitted radiation strongly depend 
on it. What is also exciting is that a black hole can emit all types of particles that 
exist in nature, independently of their spin, charge, quantum numbers or interaction 
properties, as long as their rest mass is smaller than the black hole temperature; 
therefore, long-sought but yet undiscovered particles, like the Higgs fields or super- 
symmetric particles, might indeed make their appearance in the decay spectrum of 
a black hole. The launch of the LHC, or of any other future experiment able to 
probe even higher energy regimes, deserves to be awaited with great expectations 
indeed. 
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